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significant  in  the  lock-in  amplifier  approach. 

A  DLTS  system  is  described  for  measuring  deep  levels  in  diodes  exhibiting  large 
leakage  currents.  A  capacitance  bridge  is  used  employing  the  diode  to  be  tested 
along  with  a  dummy  diode  of  similar  characteristics.  The  DLTS  spectrum  of  a  leaky 
GaAs  planar  diode  is  measured  and  compared  to  experimental  results  obtained  with 
two  standard  DLTS  systems.  It  is  shown  that  measurements  with  the  standard 
systems  are  impossible  in  certain  temperature  ranges  because  of  overloading 
problems . K^he  approach  described  here,  however,  gives  the  DLTS  spectrum  between 
77^Cand  30(TK.  _ _  - 

In  the  study  or  the  main  electron  trap  (0.83eV  from  the  conduction  band)  in 
VPE  n-GaAs,  we  have  examined  the  annealing  behavior  of  this  deep  level  and  the 
effect  of  the  implanted  species.  In  the  annealing  study,  we  find  that  the  trap 
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we  observed  that  the  trap  concentration  changes  depending  on  whether  the  implanted 
species  is  incorporated  into  Ga  or  As  substitutional  sites.  The  DLTS  results 
suggest  that  the  main  electron  trap  involves  the  Ga  vacancy. 
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caused  by  the  radiation- induced  damage  are  observed  in  the  N-implanted  samples. 
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System  effects  and  data  analysis  for  deep  level  transient  spectro¬ 
scopy  (DLTS)  have  been  examined  and  applied  to  study  the  deep  levels  in 
the  GaAs-GaP  system. 

Studies  of  typical  DLTS  systems  using  either  the  lock-in  amplifier 
or  the  dual-channel  boxcar  averager  are  presented.  The  effects  of  non¬ 
zero  gate  width  for  the  boxcar  averager,  phase  angle  adjustment  for  the 
lock-in  amplifier,  and  response  time  of  a  typical  commercial  capacitance 
meter  are  investigated.  Errors  introduced  in  the  measurements  by  these 
effects  are  calculated  for  typical  cases.  Measurements  of  the  gold  level 
in  silicon  are  presented,  along  with  calculated  corrections.  We  find  the 
correction  to  be  minimal  in  the  boxcar-averager  method,  but  significant 
in  the  lock-in  amplifier  approach. 

A  DLTS  system  is  described  for  measuring  deep  levels  in  diodes 
exhibiting  large  leakage  currents.  A  capacitance  bridge  is  used  employing 
the  diode  to  be  tested  along  with  a  dummy  diode  of  similar  characteristics. 
The  DLTS  spectrum  of  a  leaky  GaAs  planar  diode  is  measured  and  compared 
to  experimental  results  obtained  with  two  standard  DLTS  systems.  It  is 
shown  that  measurements  with  the  standard  systems  are  impossible  in 
certain  temperature  ranges  because  of  overloading  problems.  The  approach 
described  here,  however,  gives  the  DLTS  spectrum  between  77°  and  300°K. 
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In  the  study  of  the  main  electron  trap  (0.83eV  from  the  conduction 
band)  in  VPE  n-GaAs,  we  have  examined  the  annealing  behavior  of  this  deep 
level  and  the  effect  of  the  implanted  species.  In  the  annealing  study, 
we  find  that  the  trap  concentration  change  is  related  to  the  Ga  vacancy. 

In  the  implantation  study,  we  observed  that  the  trap  concentration  changes 
depending  on  whether  the  implanted  species  is  incorporated  into  Ga  or  As 
substitutional  sites.  The  DLTS  results  suggest  that  the  main  electron 
trap  involves  the  Ga  vacancy. 

Electron  traps  in  VPE  n-GaAs,  P  alloys  with  several  compositions 
have  been  examined.  It  is  observed  that  the  DLTS  spectra  for  the  as- 
grown  alloys  are  similar  to  each  other.  There  are  differences  in  DLTS 
spectra  among  the  N-doped,  N-free,  and  N- implanted  samples  for  the 
same  alloy  composition.  Deep  levels  caused  by  the  radiation- induced 
damage  are  observed  in  the  N-implanted  samples. 
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I .  INTRODUCTION 

Defects  and  deep  impurity  centers  in  semiconductor  materials 
are  sites  of  generation,  recombination,  and  trapping  of  electrons  and 
holes.  These  centers  have  important  effects  on  the  electrical  and 
optical  performance  of  semiconductor  devices.  Some  of  these  centers  are 
introduced  intentionally  during  the  device  fabrication  processes  to  pro¬ 
duce  desirable  electrical  or  optical  characteristics,  for  example,  gold 
centers  are  used  in  Si  switching  devices,  chromium  in  semi -insulating 

GaAs,  and  nitrogen  isoelectronic  traps  in  GaAs,  p  LED's.  However,  some 

*"X  x 

deep  levels  are  incorporated  unintentionally,  and  can  undesirably  lower 
the  yield  of  good  devices  and  circuits ,  such  as  nonradiative  recombination 
centers  in  GaP  LED's  and  the  radiation  damage  after  ion  implantation. 

The  generation-recombination-trapping  kinetics  associated  with 
a  single  level  in  the  bandgap  have  been  analyzed  by  Shockley  and  Read 
[1],  Hall  [2],  and  extended  to  multilevel  systems  by  Sah  and  Shockley 
T3].  Deep  levels  in  semiconductors  have  been  poorly  understood  due  to 
the  lack  of  convenient  and  effective  experimental  methods,  as  well  as 
the  theoretically  difficult  nature  of  this  topic  [4] .  A  transient  ca¬ 
pacitance  technique  [5]  has  been  developed  in  conjunction  with  thermally 
stimulated  current  (TSC)  and  photoconductivity  measurements  to  provide 
information  on  the  electrical  properties  of  deep  impurities  and  defects 
[6,7].  More  recently,  a  powerful  transient  capacitance  spectroscopic 
method,  called  Deep  Level  Transient  Spectroscopy  (DLTS),  has  been  de¬ 
veloped  by  Lang  [8],  The  DLTS  method  is  particularly  useful  for  studying 
levels  deeper  than  about  O.leV  in  the  bandgap,  whether  they  are  radiative 
or  not  f4].  DLTS  is  therefore  a  complementary  tool  to  photoluminescence 


measurements,  which  are  suitable  for  studying  only  radiative  centers  [4] . 
This  technique  has  been  applied  to  the  study  of  bulk  deep  levels  as  well 
as  to  interface  states  [8-12] . 

Ion  implantation  has  been  accepted  as  one  of  the  most  important 
techniques  for  device  fabrication,  as  well  as  for  semiconductor  research. 
Ion-implanted  semiconductor  material  contains  considerable  lattice  damage 
[13,14],  and  post-implant  annealing  is  necessary  to  obtain  high  electrical 
and  optical  activation.  Recently  laser-beam  [15-19]  and  electron-beam 
[20,21]  techniques  have  been  used  to  supplement  or  replace  conventional 
high  temperature  thermal  annealing  in  the  reordering  of  implanted  ma¬ 
terial.  However,  these  two  annealing  techniques  are  still  in  their  in¬ 
fancy.  In  this  work,  ion  implantation  is  used  as  a  powerful  tool  f . r 
its  doping  controllability,  uniformity  and  reproducibility,  and  the  con¬ 
ventional  high  temperature  thermal  anneal  is  used  because  of  its  relia¬ 
bility  and  widespread  applications. 

Although  the  epitaxial  growth  of  III-V  compound  materials  has 
made  significant  progress,  the  existence  of  discrete  deep  states  in  con¬ 
siderable  concentrations  in  these  materials  is  an  obvious  limitation. 
Considerable  efforts  are  underway  to  assess  the  deep  level  content  with 
respect  to  concentration,  electronic  properties,  and  physico-chemical 
origin.  For  our  study  of  epitaxial  GaAs  material,  we  concentrate  on  a 

dominant  electron  trap  which  is  closely  related  to  the  growth  process. 

14  -3 

This  trapping  center  is  usually  in  the  order  of  10  cm  concentration 
in  VPE  GaAs  [22,23]  and  is  not  detected  in  typical  LPE  or  MBE  material. 

For  the  GaAs^  ^P^  alloy,  there  have  been  limited  reports  on  deep  levels 
by  the  DLTS  method  [  24-28] .  In  this  work,  we  examine  the  deep  levels 


in  various  available  alloy  compositions  of  the  ternary  alloy.  This  is 
an  important  step  for  future  work  in  this  laboratory  to  understand  the 
behavior  of  deep  levels  in  variuos  alloy  composition,  with  the  eventual 
goal  of  developing  a  unified  theoretical  explanation  of  deep  level 
behavior. 

A  general  background  review  is  given  in  Chapter  2.  The  DLTS 
setup,  system  effects,  and  data  analysis  are  addressed  in  detail  in 
Chapter  3.  The  experimental  results  and  discussions  concerning  the  work 
on  GaAs  and  GaAs^^P^  are  elaborated  in  Chapter  4.  Finally  a  summary  and 
suggestions  for  future  study  are  given  in  Chapter  5. 


2 .  BACKGROUND  REVIEW 


2.1  Generation-Recombination-Trapping  Mechanism 


The  three  most  important  processes  by  which  transitions  occur 
at  deep  levels  are  thermal,  optical,  and  Auger  processes.  For  the 
thermal  process,  energy  and  momentum  are  conserved  with  participation 
of  phonons;  for  the  optical  process  energy  conservation  is  provided  by 
photons,  and  momentum  conservation  by  phonons;  in  the  Auger-impact 
process,  other  carriers  accommodate  both  energy  and  momentum  conserva¬ 
tion.  These  processes  are  shown  for  a  single  level  in  the  bandgap  in 
Fig.  2.1,  where  [6] 

E^,  Ey  ■  conduction  and  valence  band  edges,  respectively. 

ET  »  energy  position  of  the  deep  level. 

n,  p  “  electron  and  hole  concentrations,  respectively. 


=  occupied  and  empty  deep  level  concentrations, 
respectively. 

=  electron  thermal  and  optical  capture  rates,  respectively. 

*  electron  Auger  capture  rates  by  excitations  of  another 
electron  and  hole,  respectively. 

*  electron  thermal  and  optical  emission  rates  respectively. 
■  electron  impact  emission  rates  by  impact  with  another 

electron  and  hole,  respectively. 


The  parameters  c  and  e  with  subsoript  p  are  analogous  quantities  for  holes 


2.2  Junction  Capacitance  Transient  Technique 


5 


abed  abed  ababeded 

TRANSITION  RATE 


a 

electron  capture 

(Co  *•  <3  »  ejn  ♦  cgp)  npT 

b 

ELECTRON  EMISSION 

(ej,  *  «g  »  egn  *  egp )  nT 

c 

hole  capture 

(Cp  *  eg  *  cgp  *  cgn)  pnT 

d 

hole  emission 

(«p  ♦  «P  *  egp  ♦  egn)  pT 

RATE  EQ  IS 

drvr 

>  a-b-c*d 

O  T 

Fig.  2.1.  Electron  and  hole  capture  and  emission  transitions;  the  trans¬ 
ition  rate  of  the  thermal,  optical  and  Auger-impact  processes 
at  deep  level  (Sah  [6]). 
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where  V  “  °T  +  PT 

The  advantage  of  studying  deep  levels  in  a  reverse-biased  p-n 
junction  or  a  Schottky  barrier  is  the  linearization  of  rate  equations. 
Also,  data  accuracy  is  better  compared  with  low  temperature  Hall 
measurements  [29].  In  the  depletion  region. 


n  =  p  ~  0 


dnT 

dt  *  ^en  +  ep  ^  “t  +  ep  NTT 


and 


Vc>  - 


e„  +  e„ 
p  n 


NTT 


+ 

P  n 


exp  [ -(ep  +  eQ)t] 


(2.2) 


if 


VO) 


N. 


TT 


The  thermal  emission  rate  e  can  be  related  to  the  thermal 

n 

capture  rate  c  and  thermal  capture  cross  section  a  by  detailed 
n  n  ' 

balance  at  thermal  equilibrium: 


c  N  g  e"(EC_ET)/kBT  -  a  <v  >N  g  e'(EC_ET)/kBT  (2.3: 
nun  n  n  c  n 
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where  N^,  and  are  the  effective  densities  of  states  in  the  conduction 
band  extrera;  gR  and  g^  are  the  degeneracy  factors  of  the  deep  level; 
<vr>  and  are  the  root-mean-square  (rms)  thermal  velocities  of 


electrons  and  holes;  and  k  is  the  Boltzmann  constant.  Generally,  for 

an  energy  level  in  the  upper  (lower)  half  of  the  bandgap,  the  electron 

(hole)  emission  rate  is  much  higher  than  the  hole  (electron)  emission 

rate,  so  that  the  excitation  of  trapped  holes  (electrons)  can  be 

neglected.  A  deep  level  with  e  >5*e  (e  »e  )  is  classified  as  an 

n  p  p  n 

electron  (hole)  trap.  A  deep  level  with  e^e  is  a  recombination 
center. 

For  an  electron  trap,  Eq.  (2.2)  can  be  simplified  to 
n^t)  -  exp(-ent),  if  1^(0)  - 

As  shown  in  Fig.  2.2  and  Fig.  2.3,  electron  traps  which  are  majority 
carrier  traps  can  be  filled  up  initially  by  a  voltage  pulse  with  an 
adequate  duration,  which  collapses  the  depletion  region.  Hole  traps 
which  are  minority  carrier  traps  can  not  be  filled  unless  the  voltage 
pulse  is  large  enough  to  forward  bias  the  junction  and  inject  the 
minority  carriers. 

It  is  known  that  the  depletion  region  capacitance  of  an 

*|> 

abrupt  p-n  junction  or  a  Schottky  barrier  on  an  n-type  semiconductor 
is  given  by 


where  Vbi  is  the  built-in  voltage  of  the  junction  or  barrier;  VR  is  the 
reverse  bias;  €  is  the  permittivity  of  the  semiconductor  material;  q  is 
the  electronic  charge;  ND+  is  the  net  charge  of  uncompensated  donors, 
which  are  assumed  constant  in  the  depletion  region;  A  is  the  area  of  the 
junction  or  Schottky  barrier;  W  is  the  depletion  width.  An  analogous 
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Fig.  2.2.  Capacitance  transient  due  to  a  majority  carrier  trap  in  a  p"**-n 
diode.  The  insets  labeled  1  through  4  schematically  show  the 
charge  state  of  the  defect  level  and  width  of  the  space  charge 
region  (unshaded  portion)  at  various  times  before  and  during 
the  transient  (Lang  [8J). 


Decay  of  Transient 
Due  to  Thermal  Emission 
t>0 

L  P-1478 


2.3.  Capacitance  transient  due  to  a  minority  carrier  trap  in  a 
p  -n  diode  (Lang  [8]). 


expression  holds  for  an  abrupt  n  p  or  Schottky  barrier  on  p-type  semi¬ 
conductor  material. 


As  shown  in  Fig.  2.2,  before  the  trap-filling  pulse  there  is 
a  quiescent  depletion  width.  The  rapid  and  large  capacitance  change  due  to 
the  depletion  width  reduction  during  the  pulse  period  sometimes  overloads 
the  measurement  system.  Immediately  after  the  pulse,  the  carriers  in 
the  shallow  donor  states  are  thermally  excited  into  conduction  band 
and  swept  away  by  the  junction  electric  field.  The  capacitance  then 
decreases  to  the  value  contributed  by  the  net  ionized  donors  and  the 
electrons  in  the  trapping  centers.  Thermal  emission  of  the  trapped 
electrons  is  a  slow  process,  and  the  time  constant  tr  (l/eQ)  depends 
mainly  on  the  ratio  of  the  activation  energy  to  the  device  temperature. 
Therefore,  a  slow  capacitance  transient  can  be  observed  in  a  suitable 
temperature  range  and  utilized  to  extract  the  important  parameters 
characterizing  this  deep  level.  If  we  assume  the  trapping  centers  are 
neutral  when  filled  with  electrons,  then  after  the  pulse  the  total 
charged  centers  are 

N+  -  Nd+  +  Ntt  (l-e't/Ta) 

The  capacitance  transient  can  be  obtained  mathematically: 


C(t)  -  A  +NTT  (1'e 


&C  -  CO  -  C(0)  -  CC-)  riV  N  +  +  N  ]  (2.7 

D  TT 

The  exact  trap  concentration  in  the  depletion  region  can  be  obtained 
fromAc  if  N_  is  known.  If  TT  <0.1,  Eq.  (2.6)  can  be  simplified  as 
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By  using  the  junction  capacitance  transient  technique,  the  trap  concen- 
tration,  the  thermal  emission  time  constant,  and  hence  the  thermal 
activation  energy  can  be  obtained  ideally  [29].  For  this  method,  the 
existence  of  only  a  single  level  must  be  known  in  advance.  For  a  multi¬ 
level  system  it  is  not  practical  to  use  the  computer  to  extract  the 
time  constants  from  the  capacitance  transient  data.  The  DLTS  method 
invented  by  Lang  [8]  is  more  powerful  and  is  becoming  increasingly 
popular  in  semiconductor  research. 

2.3  Deep  Level  Transient  Spectroscopy 

DLTS  is  basically  a  pulsed  junction  capacitance  transient 
technique.  The  most  important  feature  is  the  capability  of  displaying 
a  spectrum  of  deep  levels  as  a  function  of  device  temperature.  The 
spectrum  indicates  the  number  and  type  (majority  or  minority  carrier)  of 
traps.  Trap  concentration,  energy  level  and  capture  cross  section  can 
also  be  obtained  to  characterize  the  deep  level. 


Consider  a  reverse  biased  p  -n  junction.  A  periodic 
capacitance  transient  is  obtained  if  the  junction  is  pulsed  periodically 
to  fill  the  traps.  Since  the  thermal  emission  rate  depends  strongly 
on  the  ratio  of  the  activation  energy  to  the  device  temperature,  the 
time  constants  for  different  energy  levels  may  vary  by  orders  of 
magnitude  at  a  given  same  temperature.  The  key  concept  of  the  DLTS 
method  is  to  preset  a  rate  window  electronically  such  that  only  the  trap 
with  an  emission  rate  within  the  rate  window  can  be  observed  by  the 
DLTS  system.  The  junction  is  thermally  scanned  so  that  each  trap  can 
be  observed  at  the  temperature  where  its  emission  rate  is  within  the 
chosen  rate  window.  Either  a  dual-channel  boxcar  averager  or  a  lock-in 
amplifier  [30-32]  can  be  used  to  set  the  rate  window  and  process  the 
periodic  capacitance  transient.  The  details  of  this  method  will  be 
discussed  in  Chapter  3. 

Figure  2.4  is  an  illustration  of  how  a  dual-channel  boxcar 
averager  is  used  to  define  the  rate  window.  The  boxcar  averager  measures 
the  difference  of  the  capacitance  transient  at  two  sampling  instants 
t^  and  t2 .  The  boxcar  output  is  thus  given  by 

S (t  )  -  C(tx)  -  C(t2)  -  AC(e"tl/T-e_t2/T)  (2.10) 

The  signal  S(t)  will  reach  a  maximum  at  the  temperature  where  the 

thermal  emission  time  constant  t  is  obtained  from  ^ ^  I  *  0. 

max  dT  T 

max 

°r  Tmax  *  ‘  As  discussed  above,  the  DLTS  peaks  for 

various  levels  appear  at  specific  temperatures,  and  a  DLTS  spectrum  is 
displayed  in  a  single  thermal  scan. 
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Fig.  2.4.  Illustration  of  how  a  dual-channel  boxcar  averager  is  used  to 
define  the  rate  window.  The  left  hand  side  shows  capacitance 
transients  at  various  temperatures,  while  the  right-hand  side 
shows  the  corresponding  DLTS  signal  resulting  from  using  the 
boxcar  averager  to  display  the  difference  between  the  capa¬ 
citances  at  times  t.  and  t-  as  a  function  of  temperature 
(Lang  [8]).  L  z 


Temperature 


As  shown  in  Fig.  2.2  and  Fig.  2.3,  the  capacitance  transients 
due  to  a  majority  carrier  trap  and  a  minority  carrier  trap  are  of 
opposite  polarity,  as  are  the  DLTS  peaks.  Therefore,  a  DLTS  spectrum 
can  indicate  clearly  the  type  of  the  trapping  centers. 

2.3.2  Trap  concentration 

The  DLTS  peak  is  proportional  to  AC,  which  is  related  to  the 
trap  concentration  as  given  by  Eq.  (2.7)  or  Eq.  (2.9). 

2.3.3  Trap  energy  level 

Equation  (2.3)  can  be  written  as 
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Where  m  is  the  conductivity  electron  effective  mass;  is  the  density 


of  states  electron  effective  mass;  h  is  Planck's  constant;  and  E^-E^, 


is  the  energy  depth  of  the  trap  below  the  conduction  band  minimum.  It 
has  been  shown  that  nonradiative  capture  and  recombination  by  multi¬ 
phonon  emission  [33]  (MPE)  via  lattice  relaxation  instead  of  cascade 
capture  [34]  is  a  commonly  occurring  nonradiative  mechanism  in  both 
GaAs  and  GaP  [33].  At  sufficiently  high  temperature,  the  MPE  process 
causes  the  capture  cross  section  to  increase  exponentially  with 


temperature.  It  is  proposed  that  [24,33] 
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here  is  obtained  from  thermal  capture  rate  measurement  as  a  function 


of  temperature.  Therefore 
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where 


AE  “  Ec  -  Ej  +  Wa 

Therefore, 
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From  the  preset  emission  rate  windows  and  the  corresponding  temperatures 

where  the  DLTS  peaks  for  this  level  occur,  an  Arrhenius  plot  of 
2 

T  VS  l/T  can  be  obtained.  The  energy  depth  aE  is  thus  obtained  from 
the  slope  of  this  plot. 

2.3.4  Capture  cross  section 

The  magnitude  of  the  carrier  capture  cross  section  depends 
strongly  on  the  capture  process,  charge  states  of  the  trapping  centers, 
temperature,  and  electric  field  in  the  region  of  interest  [35].  Capture 
cross  section  can  be  calculated  from  Eq.  (2.13)  or  can  be  measured  by 
the  DLTS  method  as  discussed  here. 

During  the  trap-filling  pulse  period,  the  parameters  in 
Eq.  (2.1)  are  of  different  values  from  those  under  a  quiescent  reverse 
bias.  For  an  dlectron  trap  in  the  majority  carrier  case. 


e  *  c  n.  «  c  n  if  E  -E_  >**  E  -E 
n  n  i  n  Li  L*r 


nT(t)  -  Ntt  (1-e  n  if  nT(0)  -  0 

where  at  the  end  of  pulse  duration  w,  the  filled  trap  concentration  is 


n^w)  -  Ntt  (1-e  cntIw^ 

n_(w)  ■  N__  if  w  »  (c  n) 
T  TT  n 


(2.15) 


By  observing  the  DLTS  peak  values  for  different  pulse  durations. 


n^,  (w)  replaces  NTT  in  Eq.  (2.9): 

nT(w) 


(2.16) 


The  thermal  capture  rate  can  be  obtained  from  Eq.  (2.15)  and  Eq.  (2.16). 

The  thermal  capture  cross  section  can  then  be  calculated  from 

c  a  a  <  v  >.  In  order  to  do  the  capture  cross  section  measurement, 
n  n  n 

the  pulse  generator  should  have  a  rise  time  as  fast  as  possible  because 
(cRn)  ^  is  in  the  order  of  nanoseconds  to  microseconds. 

2 .4  Ion  Implantation 

Ion  implantation  is  now  a  standard  technique  of  doping  semi¬ 
conductors  in  the  areas  of  field-effect  transistors,  integrated  circuits, 
microwave  devices,  and  optical  detectors.  It  has  the  advantage  of  doping 
uniformity,  controllability,  and  reproducibility  when  compared  with 
conventional  thermal  diffusions.  Furthermore,  ion  implantation  can  be 
utilized  to  introduce  special  ion  species  into  semiconductors  to  produce 
desired  properties.  For  instance,  in  studies  of  defects  or  deep  impuri¬ 
ties  in  semiconductors,  one  can  implant  certain  ion  species  to  create 
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deep  levels  and  characterize  these  levels  in  terms  of  activation  energies 
as  well  as  capture  cross  sections.  These  results  can  pave  the  way  to 
identify  the  deep  levels  which  are  important  to  the  understanding  of 
material  properties  and  device  physics. 

Ion  implantation  denotes  the  bombardment  of  a  solid  substrate 
with  ions  in  the  keV  to  MeV  energy  range.  A  typical  implanter,  MP  300 
keV  ion  implanter  from  Accelerators  Inc.,  is  shown  schematically  in 
Fig.  2.5.  It  consists  of  an  ion  source,  a  mass  separator,  and  a  target 
chamber.  Source  materials  are  introduced  into  the  ion  source  and  excited 
to  obtain  the  desired  ions.  These  ions  are  extracted  from  the  source 
and  accelerated  by  a  high  electric  field.  The  ion  beam  is  then  deflected 
by  a  magnet  and  slit  configuration  so  that  the  desired  ions  can  be 
separated  from  unwanted  species  having  different  ratios  of  ion  mass  to 
ion  charge.  The  ion  beam  is  scanned  over  the  target  area  to  obtain  a 
homogeneous  doping.  In  the  target  chamber,  the  sample  holder  is  tilted 
from  the  beam  normal  by  an  angle  of  7°  in  order  to  avoid  channelling  of 
ions  along  specific  crystal  directions. 

The  distribution  of  implanted  ions  in  the  solid  substrate 
depends  on  the  acceleration  voltage  and  the  masses  of  the  implanted  ions 
and  the  host  atoms.  Stopping  of  the  high  energy  ions  by  the  target  atoms 
occurs  by  interactions  with  core  electrons  and  by  nuclear  collision 
processes.  These  primary  collisions  cause  displacement  of  host  atoms, 
which  in  turn  are  stopped  by  cascade  collision  processes.  Hence,  a 
statistical  distribution  of  implanted  ions  occurs  within  the  solid.  For 
an  amorphous  target  this  distribution  is  peaked  around  the  mean  projected 
range  Rp  with  a  standard  deviation  ^R^.  Tables  of  and  ^R^  are 


Schematic  diagram  of  the  ion  implantati 
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available  [36],  computed  from  the  Lindhard-Scharff  -Schiott  (LSS) 
theory  T37].  The  implanted  ion  profile  can  usually  be  approximated  by 
a  Gaussian  function  given  by 

2 

N(x)  -  N  exp  [  ~(X~V  3  (2.17) 

P  2 

2AR  * 

P 

where  N  is  the  peak  concentration  obtained  at  x  *  R  : 

P  P 

N  -°-4NP 
P  ar 

p 

and  is  the  fluence  (also  called  dose)  of  the  implant  in  unit  of 
2 

ions/cm  .  Some  deviation  of  the  profile  from  a  Gaussian  shape  is  i 

usually  expected  because  of  higher  order  moments  [36]  neglected  in 
Eq.  (2.17). 

The  implanted  ions  are  introduced  into  the  crystal  by  brute 
force.  The  incorporation  of  foreign  atoms  into  the  crystal  lattice  is 
therefore  far  from  what  is  expected  at  thermal  equilibrium.  The  lattice 
damage  depends  on  the  dose,  acceleration  voltage,  relative  ion  mass  to 
target  atom  mass,  and  the  temperature  at  which  the  implant  is  performed 
[14].  The  location  of  implanted  atoms  in  the  lattice  may  be  either 
substitutional  or  interstitial.  Most  defects  caused  by  the  implantation 
are  unstable  [14,38].  Even  at  room  temperature  they  tend  to  form  more 
stable  configurations  [14, 30] .  In  order  to  restore  the  lattice  order  and 
semiconductor  properties  in  the  bombarded  region,  the  substrate  must  be 
annealed  at  elevated  temperatures.  Although  recrystallization  of 
amorphized  silicon  layers  can  be  achieved  by  annealing  at  ~500°C-600°C 
for  a  few  minutes  in  an  inert  ambient  [39],  residual  defects  are  observed 
by  the  DLTS  method  [40,41].  The  annealing  behavior  is  therefore  very 


complex.  Impurity  atoms  have  to  compete  with  lattice  atoms  and  defects 
to  obtain  a  stable  incorporation.  In  order  to  fully  recover  the 
electrical  and  optical  properties,  annealing  temperature  as  high  as 
900  ~  1000°C  are  necessary.  For  compound  semiconductors  protective 
surface  layers  or  suitable  ambients  are  necessary  to  avoid  decomposition 
of  the  material  due  to  the  incongruent  vaporization  of  compound  elements 
[42].  Encapsulations  by  Si^N^  or  Si02  layers  are  used  in  this  work. 

2.5  Encapsulation  and  Anneal 

2.5.1  Silicon  dioxide  deposition 

Shown  in  Fig.  2.6  is  the  SiC^  deposition  system,  in  which 
SiO^  films  are  deposited  by  the  oxidation  of  silane.  A  mixture  of  silane, 
nitrogen  and  oxygen,  flow  over  the  sample  mounted  on  a  heated  plate  at 
500°C.  The  flow  of  gases  is  controlled  by  needle  valves  and  monitored 
by  flowmeters.  The  approximate  flow  rates  are  ~  80  standard  cubic 
centimeters  per  minute  (SCCM)  of  nitrogen,  ~  1.1  SCCM  of  oxygen,  and 
~  28  SCCM  of  2%  silane  with  98%  nitrogen.  Films  used  in  this  work  are 
typically  1000  ~  1500  &. 

2.5.2  Silicon  Nitride  Deposition 

Figure  2.7  shows  schematically  the  RF  plasma  Si^N^  deposition 
system  [43].  The  details  of  this  system  have  been  discussed  elsewhere  [44] 
A  brief  operation  procedure  is  given  here.  After  loading  the  samples 
on  a  graphite  heater,  the  chamber  and  gas  line  are  pumped  down  to 
<  3  x  10  ^  torr.  A  preliminary  nitrogen  discharge  is  initiated  to 


remove  remaining  oxygen  from  the  system  while  the  sample  is  covered  by 
a  shutter.  After  running  the  nitrogen  discharge  for  a  few  minutes,  the 
rf  power  is  turned  off.  After  introducing  nitrogen  and  silane  to  a 


Fig.  2.7.  Schematic  diagram  of  rf  plasma  silicon  nitride  deposition 
system  in  which  gases  are  introduced  into  reaction  chamber 
separately  (Helix  [43]). 
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desired  flow  rate  monitored  by  mass  flowmeters,  the  shutter  is  removed, 
the  sample  is  quickly  heated  to  ~  300-350°C,  and  plasma  discharge  is 
then  initiated.  For  most  encapsulation  work  the  flow  rates  are  50  SCCM 
for  Nj  and  ~  11-12.6  SCCM  for  the  2%  SiH^  in  Ar.  The  films  deposited 
are  usually  ~  1000-1500  %.  It  has  been  shown  [43,44]  that  Si^N^  films 
can  be  obtained  with  stoichiometric  ratio  (Si  :  N  *  0.75),  and  both 
Rutherford  backscattering  and  Auger  electron  spectroscopy  show  that 
the  film  contains  little  oxygen.  This  is  an  important  result,  since 
oxygen-contaminated  nitrides  allow  gallium  outdiffusion  during  the 
annealing  of  GaAs  [45] . 

2.5.3  Anneal 

tl 

The  high  temperature  anneals  are  performed  in  a  12  Trans 
Temp  furnace  with  a  silica  liner.  The  ambient  during  the  anneal  is 
flowing  N2  gas,  and  the  temperature  is  monitored  with  a  Chrome 1-Alumel 
thermocouple  and  Fluke  2100  A  digital  thermometer. 


Shown  in  Fig.  3.1  is  the  block  diagram  of  a  DLTS  system  with  a 
Boonton  72B  capacitance  meter  used  early  in  this  work.  The  pulsed  bias 
is  generated  by  an  HP-6115A  precision  power  supply  and  a  Systron-Donner 
SD114A  pulse  generator.  The  voltage  pulse  is  coupled  to  the  diode  through 
a  pulse  transformer  as  shown  in  the  diagram.  A  PAR  HR-8  lock-in  amplifier 
is  used  to  process  the  capacitance  transient  and  set  the  rate  window. 

The  lock-in  output  is  plotted  against  temperature  on  an  x-y  recorder. 

A  Sulfrian  cryogenic  dewar  was  modified  for  electrical  measurement.  An 
Artronix  temperature  controller  is  used  for  temperature  selection  and 
thermal  scanning.  Both  a  Chromel-Alumel  thermocouple  and  a  calibrated 
Rosemount  temperature  sensor  are  placed  adjacent  to  the  sample  for 
accurate  temperature  measurement. 

Distortion  of  the  pulse  waveform  can  occur  due  to  the  pulse 
transformer.  This  distortion  problem  is  critical  to  the  capture  cross 
section  measurement  where  a  fast -rise  pulse  is  required.  The  voltage 
pulse  can  be  connected  to  the  LO  bias  terminal  of  the  Boonton  72B  meter 
with  a  50  0  resistor  to  replace  an  internal  shunt  capacitor  of  0.47  p,F. 
With  this  change  in  the  DLTS  set-up  as  shown  schematically  in  Fig.  3.2, 
one  can  avoid  the  pulse  distortion  considerally  without  using  the  pulse 
transformer. 

A  typical  commercial  capacitance  meter  such  as  Boonton  72B  has 
a  response  time  of  about  1  ins.  This  relatively  slow  response  time  and 
consequent  overloading  problem  of  the  meter  during  the  pulsing  period 


25 


Fig.  3.1.  Block  diagram  of  a  DLTS  system  using 
meter  and  a  pulse  transformer. 


Boonton  72B  capacitance 


Fig.  3.2.  Block  diagram  of  a  DLTS  system  using  a  Boonton  72B  capacitance 
meter. 
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make  it  necessary  to  gate  off  the  signal  during  the  initial  period  of 
the  transient  in  order  to  avoid  spurious  electrical  effects  due  to  the 
bias  pulse.  A  detailed  discussion  of  system  effects  and  data  analysis 
[46]  is  given  in  Section  3.2. 

With  the  set-up  shown  in  Fig.  3.2,  we  have  examined  the  Au- 

doped  Si  diodes  from  the  National  Bureau  of  Standards  to  calibrate  our 

DLTS  system.  In  the  fabrication  of  these  diodes  gold  was  evaporated 

on  the  back  side  of  a  Si  wafer  and  diffused  for  24  hrs.  at  825°C  to 

13  -3 

obtain  a  gold  concentration  of  approximately  10  cm  .  The  diode 

diameter  used  in  this  work  was  17  mils.  Shown  in  Fig.  3.3  are  the  DLTS 

spectra  of  the  electron  trap  of  Au  in  a  typical  p+n  junction.  By 

adjusting  the  lock-in  frequency  to  preset  the  rate  window,  we  can  obtain 

a  DLTS  spectrum  for  each  rate  window.  From  the  emission  rates  and 

corresponding  temperatures  at  the  DLTS  peaks,  the  energy  level  of  this 

trap  is  estimated  to  be  0.54  +0.01  eV  as  shown  in  Fig.  3.4.  The  capture 

-15  2 

cross  section  O'  is  estimated  to  be  2  x  10  cm  at  room  temperature, 
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and  the  trap  concentration  is  about  2  x  10  cm  .  These  data  are 
consistent  with  previous  work  [47,48]  and  our  system  is  properly  calibra¬ 
ted.  The  decrease  of  the  DLTS  peak  for  higher  frequency  will  be  fully 
discussed  in  Sec.  3.2. 

3.1.2  DLTS  set-up  with  a  specially  designed  capacitance  bridge 
In  the  junction  capacitance  transient  technique,  a  capacitance 
transient  is  obtained  due  to  the  "slow"  thermal  emission  process  of 
carriers  in  the  trapping  centers.  The  emission  is  slow  relative  to  the 
period  of  the  driving  signal  of  a  capacitance  bridge.  The  effect  of 
driving  frequency  on  the  capacitance  transient  due  to  trapping  centers 
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.  DLTS  spectra  of  the  electron  (majority  carrier)  trap  in 
the  Au  diffused  p+-n  silicon  junction  for  various  lock-in 
frequency  settings. 
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has  been  experimentally  demonstrated  [49].  A  wider  range  of  trap  energy 
and  better  S/N  ratio  can  be  obtained  from  higher  driving  frequency.  Dr. 

D.  V.  Lang  proposed  a  DLTS  system  with  a  capacitance  bridge  of  20  MHz 
[50]  instead  of  1  MHz  as  used  in  Boonton  72B.  As  shown  in  Fig.  3.5 
is  the  block  diagram  of  an  improved  DLTS  system  with  a  homemade  capaci¬ 
tance  bridge.  Tektronix  SG  503  provides  a  constant  rf  voltage  output 

of  20  MH  .  The  rf  signal  power  is  split  into  the  reference  and  test 
z 

paths  without  interaction  between  the  two.  The  reference  rf  voltage 
should  be  as  large  as  possible  to  drive  the  HP  rf  mixer  into  saturation. 

In  the  testing  path,  the  rf  voltage  across  the  test  diode  is  kept 
<  0.2V  peak  to  peak  [50]  by  a  variable  attenuator  and  the  50  ft  output 
impedance  of  the  pulse  generator.  The  diode  quiescent  capacitance  is 
nulled  out  by  the  attenuating  and  phase-shifting  circuit  in  the  dummy 
arm.  Any  impedance  change  of  the  test  diode  can  be  sensed  and  amplified 
by  the  PAR  115  low  noise  wideband  amplifier.  \This  signal  is  demodulated 
by  the  reference  signal  through  a  phase  sensitiveNJatector .  The  phase 
shifter  in  the  reference  path  is  used  to  pick  up  either  the  conductive 
or  the  capacitive  component  of  the  impedance  change  in  the  test  diode. 

A  Boonton  76-3A  standard  capacitor  across  the  diode  is  used  to  discrimi¬ 
nate  the  capacitive  mode,  by  maximizing  the  PAR  113  low  noise  pre¬ 
amplifier  output  in  response  to  a  1  pF  change  in  the  standard  capacitor. 
Although  there  are  two  other  methods  for  mode-selection  proposed  by  Dr. 

Lang,  this  is  the  most  straightforward.  The  total  system  response  time 
depends  only  on  the  PAR  113  output  time  constant,  which  is  usually  on 
the  order  of  1  usee.  The  system  recovery  time  from  overloading  is  <  1  us  as 
long  as  the  PAR  113  (which  has  a  long  recovery  time)  is  not  overloaded. 
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DLTS  system  with  a  fast  capacitance  bridge  of  lp,sec  response 
time  (Lang  [8]). 


Shown  in  Fig.  3.6  is  the  DLTS  spectrum  of  a  p  n  LPE  GaAs  diode 

obtained  from  Lang  to  calibrate  our  capacitance  bridge.  The  p+  layer  is 
18  -.3 

Ge-doped  to  2  x  10  cm  and  the  n  layer  is  unintentionally  doped  with 

16  *3 

a  net  donor  concentration  of  1.5  x  10  cm  .  There  are  two  hole 
(minority  carrier)  traps  with  energy  levels  at  0.44  +0.01  eV  and  0.76 

+0.01  eV  from  the  valence  band  edge.  The  estimated  trap  densities  are 

14  -3  14  -3 

6.0  x  10  cm  for  the  deeper  one  and  8.3  x  10  cm  for  the  shallower 

one,  as  measured  by  Lang  [8],  In  addition,  a  n+p  gold  diffused  silicon 

diode  is  used  to  calibrate  this  system.  As  shown  in  Fig.  3.7,  there 

are  two  hole  (majority  carrier)  traps  in  the  junction  region.  The 

estimated  energy  levels  in  this  case  are  0.33  +0.01  eV  and  0.54  +0.01 


eV  from  the  valence  band.  These  data  agree  very  well  with  those  observed 
by  Miller  et  al  [9] . 

3.1.3  DLTS  system  with  the  two-diode  method 

As  discussed  previously,  DLTS  measurements  require  a  capacitance 
bridge  with  high  sensitivity  and  fast  response  to  detect  the  usually 
small  capacitance  transient.  Since  leakage  currents  are  highly  tempera¬ 
ture  sensitive,  a  large  leakage  current  always  prohibits  a  bridge  circuit 
from  being  balanced  on  a  high  sensitivity  scale  over  a  wide  range  of 
device  temperature.  Naturally  it  is  desirable  to  keep  the  leakage 
currents  of  the  diodes  as  low  as  possible.  However,  applications  of  the 
DLTS  method  in  the  study  of  semiconductor  defects  often  involve  imperfect 
junctions.  For  example,  diodes  formed  by  ion  implantation  usually 
exhibit  large  leakage  currents  if  the  post-implant  anneal  conditions 
are  not  adequate  to  remove  all  radiation  damage.  As  a  consequence,  it 
is  difficult  to  apply  DLTS  to  ion  implanted  structures,  particularly  in 
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3.6.  DLTS  spectrum  of  hole  (minority  carrier)  traps  in  a  p+-n 
GaAs  junction. 
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the  study  of  partially  annealed  material.  In  this  section  we  describe  a 
DLTS  system  using  two  similar  diodes  in  the  bridge  circuit  to  offset  the 
leakage  currents.  We  examine  the  effects  of  diode  leakage  currents  on 
the  DLTS  measurements  made  using  this  approach  compared  with  two  other 
systems.  The  experimental  results  of  a  leaky  Be-implanted  GaAs  planar 
p*n  junction  show  that  this  new  version  successfully  eliminates  the 
problem  of  measuring  diodes  with  appreciable  leakage  currents  [51] . 

The  DLTS  system  used  in  this  work  is  shown  in  the  block 
diagram  of  Fig.  3.8.  Two  diodes  having  similar  C-V  and  I-V  character¬ 
istics  are  mounted  on  the  same  header  to  minimize  the  temperature  dif¬ 
ference  between  the  two  devices.  The  20  MHz  driving  signals  across 
these  two  diodes  are  set  180°  out  of  phase  by  a  phase  shifter,  and  kept 
below  0.2  V  peak-to-peak.  Hence,  the  quiescent  capacitances  of  the  two 
diodes  under  the  same  reverse  bias  are  well  balanced  at  any  temperature. 
As  usual,  we  apply  a  periodic  voltage  pulse  for  trap  filling  and  emptying 
to  detect  deep  levels  in  the  depleted  region  of  the  junction.  In  order 
to  pulse  the  test  diode  only,  the  pulse  must  pass  through  two  variable 
attenuators,  a  bandpass  filter,  and  a  power  splitter  before  reaching 
the  dummy  diode.  The  power  splitters  are  used  for  impedance  matching 
and  decoupling  the  interactions  between  the  two  output  arms. 

A  fast  gate-off  circuit  along  with  its  input  and  output  wave¬ 
forms  is  shown  in  Fig.  3.9.  To  gate-off  the  capacitance  transient  during 
the  initial  overloading  period,  the  input  signal  is  held  at  the  value 
present  just  before  the  overloading  period  and  then  is  sampled  for  the 
remaining  period.  The  gate-off  period  (holding  period)  is  controlled  by 
the  RC  time  constant  of  the  NE  555  integrated  circuit  timer.  The  holding 
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and  sampling  commands  should  occur  right  before  and  after  the  bias  pulse 
which  causes  the  overloading.  The  synchronization  output  of  a  Systron- 
Donner  114A  pulse  generator  is  set  ahead  of  the  output  pulse  to  trigger 
the  holding  command.  The  total  response  time  of  this  gate-off  circuit 
is  less  than  2  jisec  and  therefore  is  compatible  with  the  capacitance 
bridge. 

DLTS  measurements  were  performed  using  20  V  reverse  bias  and 

20  |j,sec  pulse  width  on  two  similar  implanted  GaAs  planar  p-n  junctions. 

These  two  leaky  diodes  were  taken  from  the  same  area  of  the  GaAs  wafer, 

and  therefore  were  subjected  to  the  same  fabrication  processes.  The 
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diodes  were  formed  by  implanting  250  keV  Be  (10  cm  fluence)  into 

n-type  VPE  GaAs  through  a  photoresist  mask.  The  sample  was  encapsulated 

in  Si-N.  and  annealed  for  30  minutes  at  800°C.  The  I-V  characteristic 
3  4 

of  the  test  diode  is  shown  in  Fig.  3.10.  Three  DLTS  systems  were  used 
to  make  measurements  on  this  diode.  The  first  system  is  shown  in  Fig. 
3.8.  The  second  one,  using  a  Boonton  72B  capacitance  meter,  is  shown  in 
Fig.  3.2.  The  third  system  is  shown  in  Fig.  3.5.  A  lock-in  amplifier 
with  the  fast  gate-off  circuit  is  used  on  all  three  systems  for  this 
experiment. 

The  solid  curve  in  Fig.  3.11  is  the  DLTS  spectrum  of  the  leaky 

diode  using  the  two-diode  system.  It  indicates  that  there  are  four  deep 

electron  traps  in  the  junction  region.  The  activation  energies  are 

0.43,  0.58,  0.86  and  1.26  eV.  The  activation  energies  for  these  levels 

2 

are  obtained  from  the  slopes  of  the  Arrhenius  plots  of  T  t  vs.  1/T.  The 
thermal  emission  rate  window,  l/r,  is  set  by  a  PAR  HR-8  lock-in  amplifier 
The  0.43,  0.58,  and  0.86  eV  levels  have  been  observed  previously  in 
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3.11.  DLTS  spectra  of  the  test  diode,  having  an  I-V  characteristic 

as  shown  in  Fig.  3.10.  The  solid  curve,  obtained  by  using  the 
two-diode  system,  indicates  that  there  are  four  deep  levels. 
The  dashed  curve,  obtained  by  using  the  system  of  Fig.  3.2 
shows  that  the  capacitance  meter  moves  off  scale  due  to  the 
large  diode  leakage  current. 


undamaged  VPE  material  [22].  The  dashed  curve  in  Fig.  3.11  is  the  DLTS 
spectrum  using  the  system  which  employs  a  Boonton  72B  capacitance  meter. 
This  meter  will  not  stay  balanced  on  the  1  pF  scale  and  runs  off  scale 
during  the  temperature  scan.  Although  the  capacitance  transient  is  only 
of  the  order  of  0.01  pF,  the  leakage  current  prohibits  the  capacitance 
bridge  from  remaining  balanced.  If  the  third  system  is  used,  a  spectrum 
similar  to  that  of  the  two-diode  system  is  obtained,  but  the  position 
of  each  peak  depends  on  the  temperature  at  which  the  bridge  was  initially 
balanced.  Also,  the  magnitude  of  each  peak  is  not  reproducible.  This 
phenomenon  can  be  explained  by  considering  the  rf  mixer  as  a  phase - 
sensitive  detector,  with  the  20  MHz  reference  signal  as  the  LO  input. 

The  bridge  output  is  amplified  and  connected  to  the  RF  terminal.  There¬ 
fore  the  RF  input  includes  both  the  unbalance  residue  and  the  capaci¬ 
tance  transient.  Since  the  output  of  the  mixer  is  a  nonlinear  function 
of  the  LO  and  RF  inputs,  the  transient  part  of  the  mixer  output  will  be 
proportional  to  the  capacitance  transient  change  of  the  diode  only  if 
the  bridge  offset  is  small.  If,  however,  the  offset  changes  significantly 
during  the  wide-range  temperature  scan,  the  transient  part  of  the  mixer 
output  is  no  longer  proportional  to  the  capacitance  transient.  Therefore, 
the  peak  position  and  magnitude  depend  on  the  temperature  at  which  the 
bridge  is  balanced  initially.  This  is  precisely  the  situation  when  a 
leaky  diode  is  being  tested. 

A  new  version  of  a  DLTS  system  employing  a  two-diode  method  has 
been  examined  in  these  experiments.  Comparisons  of  results  of  this  system 
with  those  of  two  standard  systems  show  that  the  two-diode  version  is  es¬ 
sential  for  accurate  measurements  on  diodes  having  large  leakage  currents. 


Deep  Level  Transient  Spectroscopy  (DLTS)  [8]  has  become  a 
wide-spread  method  for  studying  deep-level  impurities  and  defects  in 
semiconductors.  Since  several  variations  of  the  basic  system  are  in  use, 
it  is  worthwhile  to  compare  and  evaluate  system  effects  and  methods  of 
data  analysis  which  could  affect  the  results  obtained  by  different  tech¬ 
niques.  In  its  original  concept,  the  DLTS  system  was  proposed  using  a 
dual  channel  boxcar  averager  for  establishing  a  rate  window  [8].  Another 
possibility  is  to  use  a  lock-in  amplifier  in  place  of  the  boxcar.  In 
either  case,  the  most  important  problem  is  to  extract  the  DLTS  rate 
window  from  the  boxcar  or  lock-in  settings.  This  analysis  has  been  pre¬ 
sented  in  the  literature  only  for  rather  idealized  situations  [8,32]. 

If  a  typical  commercial  capacitance  meter  is  used  to  measure 
the  transient  capacitance,  the  relatively  slow  response  time  and  conse¬ 
quent  overloading  problems  of  the  meter  make  it  necessary  to  gate-off  the 
signal  during  the  initial  period  of  the  transient  in  order  to  avoid 
spurious  electrical  effects  due  to  the  bias  pulse.  For  the  boxcar 
averager  case,  it  is  possible  to  avoid  problems  arising  from  this  gate-off 
requirement  by  setting  the  first  sampling  instant  to  occur  after  the  gate- 
off  period.  Indeed,  a  separate  gating  circuit  is  not  even  necessary  with 
a  boxcar.  For  the  lock-in  amplifier  case,  however,  the  effect  of  the 
gate-off  period  on  the  determination  of  the  rate  window  must  be  consid¬ 
ered.  This  has  not  been  done  in  the  past  [32]. 

With  the  boxcar  averager  the  measurements  are  typically  made 
using  relatively  wide  gates,  rather  than  the  narrow  sampling  periods 
originally  proposed  [8].  This  is  because  the  signal-to-noise  ratio  (S/N) 


is  better  for  larger  gate  widths.  Therefore,  a  calculation  of  the  rate 
window  is  needed  for  wide  gates . 

In  this  work  we  will  examine  these  questions  in  detail.  We 
will  calculate  the  boxcar  rate  window  for  arbitrary  gate  widths  and 
compare  this  with  the  simple  formula  originally  proposed  [8].  We  will 
extend  the  original  lock-in  rate  window  calculations  [32]  to  consider 
the  effects  due  to  the  gating-off  period.  We  find  that  the  rate  window 
is  very  sensitive  to  the  phase  setting  of  the  lock-in.  These  effects  are 
illustrated  by  experimental  examples  for  the  case  of  the  Au-donor  level 
in  Si. 

The  DLTS  system  used  in  this  work  is  shown  in  the  block  diagram 
of  Fig.  3.2.  A  Boonton  72B  capacitance  meter  is  used  to  measure  the 
transient  capacitance,  and  either  the  dual-channel  boxcar  averager  or  the 
lock-in  amplifier  is  used  for  setting  the  rate  window. 

It  is  necessary  to  gate-off  the  first  1-2  msec  of  the  capaci¬ 
tance  signal  because  of  the  1  msec  response  time  of  the  capacitance  meter 
and  the  consequent  overloading  problems.  By  using  the  fast  response  gate 
off  circuit  shown  in  Fig.  3.9,  we  can  gate-off  the  signal  for  a  chosen 
period  and  thus  avoid  baseline  restoration  problems. 

DLTS  measurements  were  performed  using  the  n+p  gold  Jiffused 
junction  diode  with  8.0  V  reverse  bias  and  20  p,sec  pulse  width. 

3.2.1  Dual- channel  boxcar  averager 

In  a  system  employing  a  dual-channel  boxcar  averager,  the 
transient  signals  are  fed  into  the  integrators  with  gates  set  at  times 
t^  and  t£  after  the  bias  pulse.  For  very  narrow  gates  the  differential 
output  S(t)  is  simply  the  capacitance  at  t^  minus  the  capacitance  at  t^ : 


S  (t )  *  C(t,)  -  C(t„).  The  rate  window  t  ,  which  is  the  value  of  the 
thermal  emission  rate  at  the  maximum  of  S(t)  vs.  temperature  (T)  for  a 
particular  trap  level,  is  then  simply  given  by^ 

’  <W  /X"  W"  O'1) 

However,  it  is  known  [9]  that  S/N  is  improved  with  wider  gate  widths  (w) . 
In  fact,  S/N  ®  Jv,  as  would  be  expected  for  white  noise.  Thus,  the  ex¬ 
pression  for  rmay  must  account  for  a  finite  sampling  period. 

The  normalized  output  s(t)  for  the  case  of  arbitrary  gate 
widths  is  given  by 

s(t)  -  »  [ C (t L )  -  C(t2)]  /AC 

-  t(1  -  e“w/T)(e‘tl/T-e“t2/T)/w  (3.2) 

where  &C  *  C  (®)  -  C(0).  The  rate  window  is  determined  by  differen¬ 
tiating  s  (t )  with  respect  to  t  and  setting  the  resulting  expression  equal 

to  zero.  Therefore,  t  is  the  solution  of 

’  max 

e'(tl"t2)/T[l+(t1/T)-e'w/T(l  +  w/T+  ^/t)] 

-[l+(t2/T)-e"w/T(l  +w/T+  t2/r)]  -  0  (3.3) 

With  numerical  computation  as  shown  in  Appendix  B  the  normalized  emission 

time  constant  r  /t,  can  be  plotted  versus  normalized  gate  width 
max  1 

p  *  w/(t2»t^).  An  example  of  such  a  computation  is  shown  in  Fig.  3.12. 
Since  t2/t^  *s  normally  kept  constant  during  the  measurement,  results 
are  plotted  in  Fig.  3.12  for  typical  values  t2/t^  *  ^  anc* 

As  a  relatively  good  approximation  to  this  exact  calculation, 
the  midpoint  of  the  gate  can  be  used  in  Eq.  (1)  to  obtain 
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Fig.  3.12.  Normalized  emission  time  constant  (Tmax/t^)  plotted  vs. 

normalized  gate  width  p  for  the  boxcar  averager  case  from 
equation  (3.3).  Also  shown  (dashed  lines)  are  approximations 
from  equation  (3.4). 


(3.4) 


Tmax  "  (t2"ti)/ln^  (t2  +  w/2>/(t1  +  w/2)l 

The  values  of  t  /t,  obtained  from  this  expression  are  plotted  as  the 
max  1 

dashed  curves  in  Fig.  3.12.  It  is  clear  that  for  cases  chosen  here,  the 
error  of  this  approximation  is  less  tnan  107,.  The  effect  of  this 
approximation  is  to  shift  the  DLTS  Arrhenius  plot  by  less  than  2  K  with 
essentially  no  change  in  slope. 

3.2.2  Lock-in  amplifier 

The  signal  response,  tuning  procedure,  and  calculation  of  trap 
concentrations  from  peak  heights  of  the  lock-in  amplifier  output  have 
been  discussed  previously  [32]  for  the  ideal  case  in  which  the  system 
response  time  and  overloading  of  the  capacitance  meter  may  be  neglected. 
However,  the  apparent  transient  signal  which  is  generated  by  the  bias 
pulse  being  filtered  with  the  response  time  of  the  capacitance  meter 
may  easily  be  confused  with  the  true  transient  due  to  deep  levels.  It  is 
therefore  common  practice  to  gate-off  a  small  portion  t^  of  the  initial 
portion  of  the  transient  signal.  A  sample-and-hold  circuit  may  be  used 
for  this  purpose  as  shown  in  Fig.  3.9.  For  the  boxcar  case  such  a 
circuit  is  actually  redundant,  since  the  gates  may  be  adjusted  to  auto¬ 
matically  reject  the  initial  portion  of  the  transient.  As  shown  in 
Fig.  3.13(a),  the  initial  value  of  the  transient  after  gating-off, 
normalized  with  respect  to  the  initial  value  before  gating-off,  is 
exp(-t^/T).  Since  t  is  a  function  of  temperature,  the  peak  value  of  this 
transient  will  vary  during  the  temperature  scan.  Therefore,  an  analysis 
of  the  lock-in  amplifier  signal  must  include  the  gate-off  time  t^.  In 
addition,  the  phase  setting  of  a  lock-in  amplifier  fundamentally  affects 
the  signal.  The  PAR  HR-8  lock-in  amplifier  used  in  Fig.  3.2  measures  the 


Fig.  3.13.  Phase  adjustment  and  gate-off  effect  in  the  lock-in  amplifier 
case;  (a)  normalized  exponential  capacitance  transient  before 
and  after  gating-off;  (b)  fundamental  Fourier  component  of  the 
transient,  clSin(2nt/T0  +  J);  (c)  lock-in  amplifier  mixer 
weighting  function  for  phase  setting  3. 


amplitude  and  phase  of  the  fundamental  Fourier  component  of  the  transient 

shown  in  Fig.  3.13(b).  The  lock-in  output  is  the  integral  of  the  product 

of  the  square-wave  weighting  function  in  Fig.  3.13(c)  and  this  fundamental 

Fourier  component.  The  phase  of  the  fundamental  component  is  0;  correspond 

ing  to  the  time  t  *  T  (0/360)  where  T  is  the  lock-in  period.  Simi- 
°  0  o  o 

larly,  the  lock-in  phase  setting  9  corresponds  to  the  time  tg  *  Tq  (9/360). 

Let  us  now  calculate  the  fundamental  Fourier  component  of  the 
transient  shown  in  Fig.  3.13(a).  The  normalized  input  signal  after 
gating-off  is  exp(-td/r)  x  exp  [-(t  +  Tq/2)/t],  where  Tq  is  the  period  of 
the  lock-in  and  hence  of  the  applied  pulse.  The  Fourier  expansion  of  a 
function  f(t)  over  the  interval  -T  /2  <  t  <  Tq/2  is 

09 

f(t)  -  ao/2  +7  [an  cos  ^  +  bQ  sin  ^  ]  (3.5) 
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f(t)  sin  dt 
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For  the  transient  in  Fig.  3(a)  the  fundamental  coefficients  a^,  b^,  and 
c^  are 

"  exp(-VT>  (2£t)  c«p(VT)“i:l/c^r>2  +n2]  (3.10 

T 

WV  "  exp(-td/T)«[exp(-To/T)-l]/r02  +  "2]  (3.11 
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From  an  integral  of  the  fundamental  Fourier  component  shown  in 
Fig.  3.13(b),  described  by  Eqs.  (3,6),  (3.12),  and  (3.13),  with  the 
square  wave  in  Fig.  3.13(c),  the  lock-in  signal  is 

SCO  -  ”p<-Vt>C1-«p(-Vt)]  cos  @1  +  6). 


(3.14 


AV  +  "2 

We  may  consider  three  basic  modes  of  lock-in  DLTS  operation:  namely, 

(1)  maximum  signal,  (2)  truncated-transient  phase  reference,  or  (3)  bias- 
pulse  phase  reference.  The  first  mode  corresponds  to  adjusting  the  lock- 
in  phase  to  give  a  maximum  DLTS  peak  for  each  lock-in  frequency.  For 
td  ■  0,  the  phase  setting  is  independent  of  frequency  and  given  [32]  by 
®  *  -24.5°.  Modes  (2)  and  (3)  are  identical  for  td  ■  0  since  the  bias 
pulse  coincides  with  the  beginning  of  the  measured  transient.  For 
td  +  0,  however,  these  modes  are  quite  different. 


The  rate  window  for  the  maximum  signal  mode  may  be  most  readily 
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obtained  from  the  solution  of  dc^/dT  *  0.  Therefore,  from  Eq.  (3.12)  we 


see  that  t  for  this  mode  is  the  numerical  solution  to 
max 


exp  cTf)  cc^  +  D[(^)2  +  "2]  -  £  ■  r C(^)2  +  "23  -  o. 


(3.15) 


The  phase  corresponding  to  this  maximum  value  of  c^  vs  T  is  given  by 


’  tan‘12TT' 


(3.16) 


max 


The  truncated -transient  phase  reference  mode  corresponds  to 
tg  a  0  in  Fig.  3.13(c),  i.e.  the  lock-in  square  wave  zero  crossing  coin¬ 
cides  with  the  beginning  of  the  signal  after  the  gating-off  period.  The 
rate  window  may  be  obtained  either  from  db^/dT  »  0  (Eq.  3.11)  or  from 
dS/dr  ”  0  with  ®  “  0  (Eq.  3.14).  Choosing  the  latter  expression  for  more 
generality,  we  have 


2t 
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The  numerical  solution  of  this  expression  for  the  case  0  »  0  gives  t 


max 


for  this  second  mode  of  operation. 

The  bias-pulse  phase  reference  mode  corresponds  to  the  lock-in 
zero  crossing  being  set  at  the  end  of  the  bias  pulse,  i.e.  at  the  true 
t  *  0  point  of  the  transient  during  the  gated -off  portion  of  the  signal. 
We  may  use  Eq.  (3.17)  in  this  case  as  well,  with  9  chosen  so  that  tg  = 
-t^  as  shown  in  Fig.  3.13. 


In  this  work,  t,  was  chosen  to  be  1.6  msec  with  the  pulse 
a 


frequency  varying  from  5  Hz  to  80  Hz.  The  numerical  solutions  as  a 


] 


function  of  lock-in  frequency  Tq  for  normalized  emission  time  constant 

t  /T  are  shown  by  the  solid  curves  in  Fig.  3.14  for  the  three 
max  o 

methods.  The  dashed  curves  show  the  (frequency  independent)  solutions 
for  the  t^  *  0  case.  It  is  clear  that  the  errors  due  to  ignoring  the 
gate-off  effect  are  significant  for  the  first  two  methods,  especially 
at  high  frequency  where  t(j/TQ  is  not  negligible.  The  third  method  -- 
referencing  the  lock-in  phase  to  the  bias  pulse  —  is  clearly  preferable, 
with  the  frequency  dependence  of  Tmay  being  only  a  few  percent  even  for 
the  worst  case. 

The  magnitudes  of  the  DLTS  peaks  obtained  using  the  three 
methods  are  shown  as  a  function  of  lock-in  frequency  in  Fig.  3.15.  For 
td  ■  0  there  is  no  frequency  dependence,  but  for  a  non-negligible  gating- 
off  period  the  effect  is  substantial.  The  DLTS  peaks  become  uniformly 
smaller  for  increased  frequency  in  all  three  methods. 

The  behavior  of  these  solutions  as  a  function  of  lock-in 
phase  setting  for  a  fixed  frequency  (50  Hz)  is  shown  in  Figs.  3.16  and 
3.17.  Figure  3.16  is  a  comparison  of  the  last  two  methods:  the  trun¬ 
cated-transient  phase  reference  mode  where  the  lock-in  zero  crossing  is 
after  t^,  and  the  bias-pulse  phase  reference  mode  where  the  zero-crossing 

is  before  t,.  The  ideal  behavior  occurs  for  t ,  ■  0  where  the  two  methods 
d  d 

are  exactly  the  same.  As  we  saw  in  Fig.  2.14,  note  again  that  the  bias 
pulse  method  is  nearly  ideal  for  t^  ^  0,  whereas  the  method  of  referencing 
the  phase  to  the  truncated  transient  may  introduce  substantial  error. 

Since  the  magnitude  of  this  error  is  frequency  dependent,  the  net  effect 
is  to  skew  the  data  on  an  Arrhenius  plot  leading  to  incorrect  activation 
energies,  as  we  shall  see  in  the  experimental  example. 


Normalized  Thermal  Emission  Time  Constant 


Fig.  3.14.  Normalized  emission  time  constant  (Tmax/T0)  vs.  pulse  frequency 
f  for  the  lock-in  amplifier  case.  The  maximum  output  method 
refers  to  equation  (3.15);  the  truncated  transient  phase  refer¬ 
ence  method  refers  to  equation  (3.17)  with  tg  -  0;  and  the  bias 
pulse  phase  reference  method  refers  to  equation  (3.17)  with 
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Fig.  3.15.  Normalized  output  magnitude  (equation  3.14)  as  a  function  of 

frequency  f  for  the  lock-in  amplifier  case.  The  three  methods 
are  the  same  as  in  Fig.  3.14.  The  dashed  line  illustrates 
neglect  of  the  gate-off  effect. 


Normalized  Thermal  Emission  Time  Constant  (Tmax/*i; 


Fig.  3.16.  Normalized  solution  of  equation  (3.17)  for  t^  ■  0  and  t<j  *  1.6msec 
as  a  function  of  the  lock-in  phase  setting  referred  to  the 
truncated  transient  and  to  the  bias-pulse  methods  at  f  *  50Hz. 


Even  if  there  is  little  frequency  dispersion  in  the  rate 

window,  we  see  in  Fig.  3.16  that  the  proper  value  of  is  quite 

sensitive  to  the  phase  setting.  The  signal  magnitude,  on  the  other  hand, 

is  relatively  insensitive  to  the  phase  setting,  as  can  be  seen  in  Fig. 

3.17.  Thus  it  is  important  that  the  phase  be  set  by  synchronizing  the 

lock-in  waveforms  with  the  signal  or  with  the  bias  pulses  and  not  by 

merely  maximizing  the  lock-in  output  signal.  A  10%  range  in  Tmqv  will 

shift  a  typical  DLTS  peak  by  less  than  about  2  K.  If  we  can  consider 

this  an  acceptable  error,  then  according  to  Fig.  3.16  we  must  set  9  to 

within  +  6°  of  the  desired  mode  of  operation  to  be  within  this  10%  range. 

More  demanding  limits  on  set  correspondingly  narrower  limits  on  A0  . 

Schott  et.  al.  [32]  have  considered  in  some  detail  a  typical 

phase  setting  routine  for  the  maximum  signal  mode  in  the  td  *  0  case. 

Note,  as  shown  in  Fig.  3.17,  however,  that  the  required  phase  setting 

^max  ^°r  t^s  mo<*e  *n  the  ?  0  case  is  frequency  dependent.  Although 

the  shift  from  0  to  50  Hz  in  the  t,  ■  1.6  msec  case  is  only  5.4°  and  thus 

a 

within  our  acceptable  limits  of  error,  one  should  be  aware  that  in 
general  the  phase  should  be  reset  for  each  frequency  when  using  this 
mode.  Similarly,  the  phase  must  be  reset  for  each  frequency  to  maintain 
the  truncated-transient  reference  mode.  On  the  other  hand,  the  bias- 
pulse  reference  mode  need  be  adjusted  only  once  for  all  frequencies. 

This,  coupled  with  the  very  small  frequency  dispersion  of  in  this 

mode,  makes  it  the  preferred  method  of  lock-in  DLTS  operation. 

3.2.3  Experimental  example  and  discussion 

In  order  to  illustrate  the  various  modes  of  DLTS  operation 


which  we  have  just  discussed,  we  measured  the  hole-emission  spectrum  of 


Normalized  Output  Magnitude 
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the  gold  donor  level  in  a  silicon  n+p  diode  using  the  apparatus  in 

Fig-  3.2.  Arrhenius  plots  were  constructed  for  data  taken  in  the  various 

operational  modes  according  to  the  rate  window  calculation  above  using 

both  the  t,  ■  0  and  t.  ■  1.6  msec  relationships  between  t  and  T  . 

d  d  max  o 

For  this  experiment,  we  set  p  *  0.1  and  *  5  for  the  boxcar  averager 

and  performed  measurements  with  f  ■  5,  20,  50,  and  80  Hz  for  the  lock-in 

amplifier.  The  gate-off  time  was  t^  *  1.6  msec.  The  relationship 

2 

and  activation  energies  ^E  (with  T  correction)  [9]  calculated  for  all 
cases  are  shown  in  Table  3.1. 

For  the  boxcar  averager,  two  observations  can  be  made  from 
these  results : 

(1)  For  the  small  setting  of  p  and  fixed  value  of  use<*  here> 

the  error  in  Tjnax  introduced  by  the  midpoint  approximation 

method  is  negligible.  The  activation  energy  ^E,  which  can  be 

2 

calculated  from  the  Arrhenius  plot  of  T  t  vs  1/T,  is  also 
obtained  with  negligible  error. 

(2)  For  large  settings  of  p  and  fixed  tg/t^,  the  error  in  the 

approximate  TmflV  is  always  less  than  10%  (Fig.  3.12).  For 

2 

particular  values  of  p  and  tj/t^,  ^  we  Plot  T  t  bs  1/T  with 

and  without  the  midpoint  approximation,  two  parallel  lines 

result.  The  activation  energy  (proportional  to  the  slope)  is 

found  to  be  the  same,  although  the  t  values  are  not 

max 

strictly  correct  in  the  approximate  method.  As  we  noted 
earlier,  however,  the  actual  shift  in  temperature  of  the 
Arrhenius  plot  is  typically  less  than  2  K. 
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For  the  lock-in  amplifier  case,  several  observations  can  be 


made : 


(1)  A  large  error  in  the  r  values  due  to  the  neglect  of  the 

max 

-t  /t 

e  d  is  expected  from  Fig.  3.14  for  the  maximum  output  and 

2 

truncated-transient  reference  cases.  The  plot  of  T  t  vs  1000/T 

for  the  maximum  output  method  (Eq.  3.15)  are  shown  in  Fig.  3.18 

for  both  the  t ,  ■  0  and  t,  ■  1.6  msec  values  of  t  .  The 
d  d  max 

data  points  plotted  from  uncorrected  t  values  (t,  *  0) 

max  v  d  ' 

deviate  somewhat  from  a  straight  line.  The  comparison  of 
corrected  and  uncorrected  data  shows  that  in  this  case  it  is 
absolutely  necessary  to  consider  the  gate-off  effect  in  order 
to  obtain  the  true  activation  energy  and  emission  time  constant. 

(2)  For  the  bias -pulse  phase  reference  method,  it  is  expected  that 
the  error  introduced  by  using  the  uncorrected  thermal  emission 
time  constants  in  the  activation  energy  calculation  should  be 
small.  As  shown  in  Table  3.1,  this  discrepancy  is  actually 
within  the  experimental  error. 

(3)  The  peak  value  of  the  lock-in  amplifier  output  is  expected  to 
decrease  with  increasing  frequency  (Fig.  3.15).  Such  a  fre¬ 
quency  variation  is  found  experimentally.  If  the  gate-off 
effect  were  truly  negligible,  the  peak  value  should  instead  be 
independent  of  pulse  frequency.  This  illustrates  that  the 
output  magnitude  should  be  corrected  in  the  calculation  of  trap 
concentration,  especially  for  higher  frequencies. 

(4)  The  phase  adjustment  is  critical  for  the  lock-in  amplifier  as 


discussed  in  connection  with  Fig.  3.16.  The  bias -pulse 


2 

•  3.18.  Arrhenius  plot  of  T  *  may  vs.  1000/T  obtained  for  the  gold  donor 
level  in  an  n+p  silicon  junction  with  and  without  the  correction 
for  maximum- output  method. 
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reference  method  involves  much  simpler  phase  settings  than 
either  of  the  other  two  methods  and  hence  introduces  less 
error.  Since  this  method  also  needs  practically  no  correction 
to  t  ,  it  is  clearly  the  method  of  preference  for  the  lock-in 
methods . 

We  see  also  that  the  boxcar  and  lock-in  methods  of  DLTS  are 
quite  comparable,  both  in  complexity  and  in  accuracy.  The  boxcar  method 
is  not  subject  to  spurious  recovery  effects  in  the  initial  portion  of  the 
transient,  but  it  is  comewhat  complex  to  change  the  gates  for  different 
rate  windows.  Using  the  lock-in  method  with  bias-pulse  phase  reference, 
it  is  relatively  simple  to  change  rate  windows;  however,  the  phase  must 
be  accurately  adjusted  and  the  initial  recovery  problem  must  be  dealt 
with.  Both  methods  yield  comparable  S/N  if  the  boxcar  is  used  with  wide 
gates.  Thus  neither  stands  out  above  the  other  as  the  "best"  method. 

In  the  end,  the  choice  of  a  rate  window  system  must  be  based  primarily 
on  personal  preference  and  on  the  best  utilization  of  existing  equipment. 

3.2.4  Summary 

The  DLTS  systems  based  on  either  a  dual-channel  boxcar  integra¬ 
tor  or  a  lock-in  amplifier  are  analyzed  for  realistic  experimental  con¬ 
ditions.  We  present  expressions  which  may  be  numerically  solved  to 
yield  the  DLTS  rate  window  for  cases  of  arbitrarily  wide  boxcar  gates 
and  for  gating-off  of  the  initial  portion  of  the  transient  in  the  lock-in 
case.  Solutions  for  typical  experimental  parameters  show  that  the  boxcar 
rate  window  may  be  adequately  approximated  by  Lang's  original  expression 
applied  to  the  midpoints  of  the  gates.  We  find  that  the  lock-in  DLTS 
systems  are  in  general  critically  sensitive  to  the  phase  setting  and 


and  analyze  in  detail  three  possible  modes  of  operation:  (1)  maximum 

signal,  (2)  phase  referenced  to  the  beginning  of  the  truncated  transient, 

and  (3)  phase  referenced  to  the  bias  pulse.  The  first  two  methods  show 

considerable  shifts  in  the  rate  window  relationship  t  /T  as  a  function 

max  o 

of  lock-in  frequency  when  the  initial  portion  of  the  transient  is  gated- 
off  to  avoid  overload  and  recovery  time  effects.  The  third  method  is 
much  less  frequency  sensitive  and  thus  is  clearly  the  best  choice  for 
those  who  prefer  lock-in  operation.  These  concepts  are  illustrated  by 
DLTS  measurements  of  the  hole-emission  spectra  of  the  gold  donor  level  in 
silicon  n+p  diodes.  The  data  show  that  a  naive  analysis  of  the  first  two 
modes  of  lock-in  operation  yields  activation  energies  which  are  seriously 
in  error.  With  proper  analysis,  however,  all  DLTS  methods  produce 
identical  results. 
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4.  STUDIES  OF  ELECTRON  TRAPS  IN  THE  GaAs-GaP  SYSTEM 
BY  DEEP  LEVEL  TRANSIENT  SPECTROSCOPY 

4. 1  Studies  of  Electron  Traps  in  GaAs 

It  has  been  reported  [22]  that  there  are  six  electron  traps  in 

device  quality  VPE  n-GaAs.  These  are  discussed  in  two  review  papers  [52, 

53]  dedicated  to  the  signatures  of  deep  levels  in  epitaxial  GaAs.  An 

electron  trap  with  thermal  activation  energy  of  0.83eV  from  the  conduction 

band  is  the  dominant  one  and  is  consistently  reported  in  DLTS  work  [23,55] 

It  has  been  conjectured  that  this  dominant  level  is  related  to  oxygen 

impurity  atoms  on  arsenic  sites,  [56-58]  or  that  it  is  related  to  the  Ga 

vacancy  [23].  Although  there  have  been  numerous  studies  of  this  level,  a 

positive  identification  of  the  responsible  defect  has  been  lacking.  In 

this  study  we  report  detailed  experimental  results  which  will  contribute 

significantly  to  the  understanding  and  identification  of  this  deep  level. 

4.1.1  Material  preparation  and  device  fabrication 

All  of  the  VPE  n-GaAs  (100)  material  used  in  this  work  was  pro¬ 
vided  by  the  Monsanto  Company  (St.  Louis).  The  epitaxial  layer  thickness 
is  about  10p.m.  Shown  in  Fig.  4-1  is  the  net  donor  concentration  profile 
of  the  undoped  n  epitaxial  layer,  obtained  from  C-V  measurement.  The 

depth  is  calculated  from  Eq.  (2.5)  and  the  capacitance  is  obtained  from 

-3  2 

an  Al-gate  Schottky  barrier  with  an  area  of  1.30x10  cm  (16mils  dia.), 

15  -3 

The  concentration  is  assumed  to  be  a  constant  value  of  1.2x10  cm  in 
this  study.  Two  wafers  from  the  same  growth  run  (Monsanto  1-1138)  are 
used  in  this  study,  and  we  find  that  these  two  wafers  are  essentially 
identical  in  terms  of  net  donor  concentration  and  the  main  trap  concen¬ 


tration. 


Liquid  phase  epitaxial  (LPE)  n-GaAs  material  was  provided  fay 


Dr.  R.  Milano  and  Professor  G.  E.  Stillman  of  the  University  of  Illinois 

at  Ur bana- Champaign.  The  LPE  epitaxial  layer,  on  an  n+  GaAs  substrate, 
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had  a  net  donor  concentration  of  about  2x10  cm  .  The  molecular  beam 

epitaxy  (MBE)  n-GaAs  used  in  this  study  was  grown  by  Professor  H.  Morko? 

at  the  Coordinated  Science  Laboratory  in  this  university  and  had  a  net 
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donor  concentration  of  1x10  cm  . 

An  A1  layer  of  1000  A  ~  2000  A  thickness  is  evaporated  to  make 
Schottky  barriers  on  epitaxial  n-GaAs.  Ag-Sn  eutectic  alloy  is  evaporated 
on  the  n+  substrate  and  sintered  at  ~  350°C  for  30  sec  for  ohmic  contact. 
Two  Schottky  barriers  are  cleaved  from  the  same  sample  and  are  mounted  on 
the  same  TO-18  or  T0-5  header  in  the  DLTS  set-up  with  the  two-diode  method 
described  previously. 

The  I-V  characteristics  of  a  typical  Schottky  barrier  on  VPE 
n-GaAs  is  shown  in  Fig.  4.2.  The  forward  I-V  characteristic  is  shown 
with  scales  of  0.5  V/div  horizontally  and  0.1  mA/div  vertically.  The  re¬ 
verse  I-V  characteristic  has  scales  of  lOV/div  horizontally  and  5pA/div 
vertically. 

4.1.2  Electron  Traps  in  VPE  n-GaAs 

As  mentioned  in  Chapter  2,  either  a  dual-channel  boxcar  averager 
or  a  lock-in  amplifier  can  be  used  in  the  DLTS  system  to  preset  the  rate 
window.  For  signals  much  less  than  lOOraV,  the  lock-in  amplifier  is  better 
than  the  boxcar  averager  in  terms  of  S/N  ratio  due  to  the  maximum  sensi¬ 
tivity  of  lOOmV  in  the  boxcar  averager.  The  gate-off  effect  discussed  in 
Chapter  3  is  negligible  in  this  case  because  of  the  fast  response  time  of 
1  ~  2nsec  in  our  DLTS  system.  Therefore,  a  lock  in  amplifier  is  preferred 


in  this  study  and  the  bias-pulse  phase  reference  mode  (which  is  the  same 
as  the  truncated  transient  phase  mode)  is  used.  The  thermal  emission 
time  constant  and  trap  concentration  calculations  are  therefore  based  on 
this  mode.  The  lock-in  amplifier  phase  is  adjusted  to  show  a  positive 
DLTS  signal  for  a  majority  carrier  trap  and  a  negative  signal  for  a  min¬ 
ority  carrier  trap. 

Shown  in  Fig.  4.3  is  the  DLTS  spectrum  of  electron  traps  in  VPE 
n-GaAs.  The  DLTS  signal  is  expressed  in  terms  of  trap  concentration  cal¬ 
culated  from  Eq.  (2.16)  by  assuming  that  each  trap  is  uniformly  distributed 
throughout  the  observed  layer.  Here  the  symbol  V  designates  VPE  material, 
and  similarly  L  and  M  stand  for  LPE  and  MBE  later  in  this  chapter.  The 
electron  trap  V7  is  dominant  in  VPE  n-GaAs.  The  other  traps  are  present 
in  only  very  small  concentrations  and  appear  only  accidently.  Therefore, 
it  is  hard  to  conduct  a  systematic  experiment  to  understand  the  physico¬ 
chemical  origins  of  traps  other  than  V7. 

The  thermal  activation  energy  for  each  trap  is  obtained  from 
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the  slope  of  the  Arrhenius  plot  of  T  t  vs  1000/T  shown  in  Fig.  4.4.  The 
leas t-square>—  fit  program  to  calculate  activation  energy  for  either  a 
boxcar  or  a  lock-in  case  is  presented  in  Appendix  C. 

The  concentration  profile  of  the  main  electron  trap  V7  is  shown 
in  Fig.  4.5,  and  the  depth  is  obtained  from  Eq.  (2.5).  The  concentration 
in  the  observed  layer  obtained  from  two  pulses  of  different  amplitude  is 
calculated  according  to  Eq.  (A7)  in  the  appendix  of  Ref.  50. 


This  equation  can  be  simplified  to 


Fig.  4.3  DLTS  spectrum  of  electron  traps  in  VPE  n- 


if  the  net  donor  concentration  is  assumed  constant  and  the  trap  con¬ 
centration  in  the  small  layer  is  taken  as  constant.  The  symbol  6  denotes 
a  small  variation  or  change  in  the  following  parameter.  In  the  calcu¬ 
lation  the  built-in  voltage  V.  .  is  assumed  to  be  0.65  V,  calculated  from 

bi 

the  relation  [59] 

Vbi  '  VB  -  (EC  -  EFs)/e  (4‘2> 

where  VQ  is  the  barrier  height  and  Epg  is  the  Fermi- level  in  semiconductor 

As  shown  in  Fig.  4.5,  the  trap  concentration  can  be  assumed  to  be  a  con- 
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stant  value  of  1.4x10  cm  . 

The  capture  cross  section  for  each  trap  at  the  observed  temp¬ 
erature  can  be  calculated  from  the  thermal  emission  rate  measurement  by 

Eq.  (2.13).  The  capture  cross  section  for  each  trap  calculated  from  the 

* 

emission  rate  measurement  is  listed  in  Table  4.1  by  assuming  m  =  0.065  m£ 

m^n  *  0.068  me,  and  g  »  1  where  mg  is  the  electron  mass.  The  capture 
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cross  section  cra  in  Eq.  (2.13)  of  1.3x10  cm  for  trap  V7  is  in  good 
agreement  with  that  reported  by  Mircea  [22] . 

The  capture  cross  section  can  also  be  obtained  from  the  thermal 
capture  rate  measurement  as  discussed  in  detail  in  Sec.  2.3.4.  Shown  in 
Fig.  4.6  is  a  typical  example  taken  from  the  VPE  n-GaAs  sample.  The 
capture  rate  is  determined  from  the  DLTS  peak  heights  as  a  function  of  the 
majority  carrier  pulse  duration.  Here,  a  Systron- Donne r  SD-110C  pulse 
generator  with  pulse  transition  time  of  5nsec  i«  used.  A  pulse  width  of 
lpsec  is  long  enough  to  fill  up  this  trap.  In  order  to  be  sure  that  no 
competition  of  carrier  capture  from  other  traps  takes  place,  a  pulse  of 
20psec  duration  is  always  used  in  the  emission  rate  measurements  in  this 
study.  From  Eq.  (2.15)  we  estimate  that  the  carrier  capture  time  con- 
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Fig.  4.6  Typical  example  of  the  dependence  of  trap  peak  heights 
on  the  duration  of  the  majority  carrier  pulse. 


stant  for  V7  is  about  40nsec.  Therefore,  the  capture  cross  section  from 

the  capture  rate  measurement  at  the  observed  temperature  (378°K)  is  esti- 
•  16  2 

mated  to  be  4x10  cm  .  This  value  is  in  good  agreement  with  that 
reported  by  Henry  and  Lang  [33,58],  They  observed  an  exponential  temp¬ 
erature  dependence  of  the  capture  cross  section  and  estimated  the  capture 
activation  energy  to  be  ~  0.08eV.  Henry  and  Lang  explained  the  strong 
temperature  dependence  by  multiphonon  emission  (MPE)  in  which  the  energy 
is  conserved  via  lattice  relaxation.  The  capture  cross  section  on  at  the 
observed  temperature  is  estimated  from  the  emission  rate  measurement  to 
be  1.1x10  cm  from  Eq.  (2.12)  for  a  =  1.3x10  cm'1'.  Hence,  there  is 

GO 

a  large  difference  between  the  capture  cross  sections  obtained  from 

capture  rate  and  emission  rate  measurements.  Mircea  et  al.  [60]  tried 

to  explain  this  discrepancy  by  artifically  introducing  a  temperature 

coefficient  a  in  addition  to  W  such  that 

a 

AE-E  -  E  +  W  -  QfT  (4.3) 

c  i  a 

This  argument  is  not  convincing,  and  has  been  challenged  by  Majerfeld 
et  al  [61],  They  proposed  that  the  transition  occurs  between  the 
valence  band  the  the  L  minimum  instead  of  the  T  minimum.  It  has  been 
reported  [62]  that  photocapacitance  measurements  on  GaAs  under  hydro¬ 
static  pressure  suggest  the  dominance  of  optical  transitions  from  trap 
levels  to  higher  lying  L  and  X  minima  rather  than  the  T  minimum.  It  has 
also  been  reported  [63]  that  the  hydrostatic  pressure  dependence  of  AE 
for  the  0.83eV  trap  in  GaAs  as  derived  from  thermal  emission  rate  measure 
ments  is  3.8  +0.3  meV/kbar.  This  value  together  with  the  value  of 
1.7+1  meV/kbar  obtained  for  optical  transitions  from  the  trap  level  to 
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the  valence  band  [64]  gives  a  pressure  dependence  of  5.5  +  1.3  meV/kbar 
from  the  valence  band  to  the  final  transition  state.  This  suggests  that 
the  final  transition  state  is  in  the  L  minimum  rather  than  the  T  minimum. 
Majerfeld  et  al  [61]  estimated  that  the  energy  depth  from  trap  level  to 
the  T  minimum  is  about  0.54eV,  using  conventional  C-V  measurements.  It 
has  been  shown  that  the  energy  separation  between  the  L  and  T  minima  in 
GaAs  is  0.29  HhO.OleV  at  300°K  [65].  The  thermal  activation  energy  of 
0.83eV  observed  in  DLTS  emission  rate  measurements  suggests  that  the  tran¬ 
sition  is  from  the  trap  level  to  the  L  minimum.  If  this  is  true,  then 

the  effective  mass  values  m  *  0.11  m  and  m,  “  0.55  m  [33]  for  the  L 

e  dn  e 

minimum  should  be  used.  The  spin  degeneracy  g^  can  be  taken  as  1/2  [66]. 
Therefore  Eq.  (2.13)  becomes 

(m  lm  )3/2 

en  -  3.26X1021  x  W-<EC  "  V/kBTl 

(m  /mg) 

-  2.40xl015a  T2exp (-0.83/k_T) 
n  a 

From  the  observed  emission  rate  and  the  corresponding  temperature  T, 
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cm  is  estimated  to  be  1.9x10  cm  .  This  value  is  in  good  agreement  with 
that  extrapolated  at  T  *  «•  by  Henry  and  Lang  [33].  Here  we  do  not  use 
the  correcting  factor  W  for  the  capture  cross  section.  The  capture  cross 
section  is  about  the  same  as  am  obtained  by  Henry  and  Lang  [33] .  There¬ 
fore,  the  capture  cross  section  does  not  have  a  temperature  dependence. 

T>.-  apparent  temperature  dependence  of  the  capture  cross  section  obtained 
by  Henry  and  Lang  is  due  to  their  assumption  that  the  available  free 
carriers  during  capture  are  from  the  r  minimum  and  with  a  value  of  net 
doping  for  a  small  trap  concentration.  Majerfeld  et  al  [61]  proposed  that 


the  capture  and  emission  processes  are  via  the  L  minimum.  The  available 
free  carrier  concentration  in  the  L  minimum  relative  to  that  in  T  has  a 


temperature  dependence 

3/2 

nL/nT  =  (mdn/mdn)  expC-AE^/kgT) 

and 

nL  +  np  =  n 

where  the  superscripts  L  and  T  stand  for  the  L  and  F  minima  respectively, 
and  AEjj  is  the  energy  difference  between  L  and  T  minima.  Although  the 
assumption  that  the  carrier  capture  is  also  via  the  L  minimum  is  obviously 
not  the  first-order  effect,  the  growing  evidence  of  this  assumption  has 
been  reported  for  other  compound  alloys  [67,68].  Whether  the  carrier 
capture  is  from  the  F  minimum  by  the  MPE  process  or  via  the  L  minimum  still 
needs  further  study. 

4.1.3  Annealing  behavior  of  electron  trap  V7  in  epitaxial  n-GaAs 
It  has  been  reported  that  the  V7  trap  concentration  decreases 
after  heat  treatment  at  ~  750° C  for  about  30  minutes  in  an  ambient  with¬ 
out  encapsulation  [22].  Mircea  et  al.  attributed  this  phenomenon  to  out- 
diffusion  of  the  defect  responsible  for  deep  level,  and  estimated  the  out- 
diffusion  activation  energy  and  outdiffusion  coefficient  as  a  function  of 
temperature.  Shown  in  Fig.  4.7  is  the  DLTS  spectrum  of  electron  traps  in 
the  unencapsulated  VPE  sample  after  heat  treatment  in  ultra  high  purity 
N2  flow  ambient  at  750°C  for  30  minutes.  Indeed  the  trap  concentration  of 
V7  does  decrease.  The  isothermal  and  isochronal  annealing  results  averaged 
over  4-8  samples  are  summarized  in  Table  4.2.  There  is  no  significant 
change  for  samples  annealed  at  temperatures  lower  than  700°C.  We  can 
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pectrum  of  electron  traps  in  unencapsulated  VPE  n-GaAs  sample  annealed 
®C  for  1/2  hr  in  N„  ambient. 
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Table  4.2.  Average  Trap  Concentrations  in  VPE  n-GaAs  for  Various 
Annealing  Conditions. 


Encapsulation 


X 


X 


X 


X 


Anneal 

Ambient 


As  overpressure 


X 

X 

1.51xl014 

700 

10 

14 

1.26x10 

750 

10 

1.07xl014 

750 

30 

5.80xl012 

700 

30 

14 

1.48x10 

800 

30 

4.01xl014 

800 

60 

5. 40x10 14 

700 

30 

1.52xl014 

800 

30 

3.34X1014 

800 

60 

2.91xl014 

750 

30 

1.41xl014 

successfully  fabricate  Schottky  barriers  on  the  samples  annealed  at  tem¬ 
perature  below  800°C,  and  obtain  reasonable  leakage  current,  forward 
and  reverse  breakdown  voltages.  The  DLTS  spectrum  shown  in  Fig.  4.8  is 
for  a  sample  annealed  at  750°C  for  1/2  h  with  As  overpressure.  The 
samples  were  sealed  in  a  quartz  ampule  at  4*10  ^  torr  with  enough  As 
powder  to  ensure  an  As  overpressure  during  the  anneal.  The  result  shows 
that  the  trap  concentration  does  not  change  significantly.  This  suggests 
that  the  annealing  behavior  cannot  be  explained  simply  by  outdif fusion 
of  a  defect.  It  also  seems  peculiar  that  the  heat  treatment  in  the 
studies  by  Mircea  et  al.  [22]  was  not  done  with  any  dielectric  encapsu¬ 
lation.  It  is  well  known  that  compound  semiconductors  tend  to  decompose 
due  to  the  incongruent  vaporization  of  compound  elements  at  temperatures 
above  600° C  [42],  We  therefore  examined  the  annealing  behavior  of  enca¬ 
psulated  samples  to  study  these  effects. 

Shown  in  Fig.  4.9  and  Fig.  4.10  are  the  DLTS  spectra  for  the 
samples  annealed  at  800° C  for  1/2  h  with  high  quality  SigN^  and  Si02 
encapsulations,  respectively.  It  is  very  surprising  that  the  trap  con¬ 
centrations  in  both  cases  do  increase  instead  of  decrease.  Shown  in 
Table  4.2  are  the  average  trap  concentrations  after  isothermal  and 
isochronal  annealing.  It  is  obvious  that  the  annealing  behaviors  are 
opposite  for  samples  with  and  without  encapsulation.  Therefore,  it  is 
reasonable  to  explain  the  annealing  behavior  by  assuming  that  the  electron 
trap  of  0.83  eV  is  related  to  the  generation  of  Ga  vacancies  rather  than 
a  simple  outdiffusion  mechanism.  The  trap  concentration  is  apparently 
related  to  the  Ga  vacancy  concetration.  For  the  annealed  unencapsulated 
samples  the  higher  (longer)  the  anneal  temperature  (time),  the  more  As 
vacancies  are  created.  Therefore,  the  effective  Ga  vacancy  concentration 
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spectrum  of  electron  traps  in  unencapsulated  VPE  n-GaAs  sample  annealed 
>0°C  for  1/2  hr  with  As  overpressure. 
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decreases,  as  does  the  trap  concentration.  For  encapsulated  samples.  As 
outdiffusion  is  suppressed  as  the  equilibrium  vapor  pressure  is  built 
up  between  the  dielectric  layer  and  the  GaAs.  It  is  also  known  that  Ga 
tends  to  diffuse  into  the  dielectric  layer,  particularly  SiO^.  There¬ 
fore,  the  Ga  vacancy  density  will  increase  for  higher  temperatures  or 
longer  anneal  periods.  This  is  confirmed  by  the  example  of  As  over¬ 
pressure.  With  As  overpressure  the  As  outdiffusion  is  prevented  and  the 
Ga  outdiffusion  at  750°C  or  below  is  not  significant.  Therefore  no 
significant  vacancy  generation  effect  is  observed.  For  higher  temper¬ 
ature  As  overpressure  experiments,  we  cannot  fabricate  good  Schottky 
barriers  due  to  the  resulting  surface  deterioration.  This  prevents 
possible  observation  of  increases  in  trap  concentration  at  higher  an¬ 
neal  temperatures. 

Although  Ga  outdiffusion  into  the  SiC^  layers  is  much  more  likely 
than  that  into  Si^N^  layers  for  high  temperature  anneal  [45] ,  appreciable 
As  outdiffusion  has  also  been  observed  [45].  Therefore,  for  samples 
encapsulated  with  SiO^  and  annealed  at  800°C  for  either  1/2  h  or  1  h, 
the  trap  concentration  really  depends  on  the  effective  concentration 
of  the  Ga  vacancy.  This  can  probably  explain  a  larger  trap  concentration 
observed  in  the  Si^N^  encapsulated  sample  compared  with  the  SiC^ 
encapsulated  one  at  the  same  annealing  condition. 

It  has  been  reported  that  the  0.83  eV  trap  does  not  exist  in 
typical  LFE  [53,  69]  or  MBE  [70]  material.  In  fact,  LFE  GaAs  material 
is  typically  free  of  electron  traps  [53].  Since  LPE  GaAs  growth 
employs  a  Ga-rich  environment  [71] ,  traps  involving  Ga  vacancies  should 
be  suppressed.  As  shown  in  Fig.  4.11,  high-quality  LPE  material  is 
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Fig.  4.11  DLTS  spectrum  of  electron  traps  in  LPE  n-GaAs. 


observed  to  be  quite  free  of  electron  traps.  The  VPE  GaAs  growth  process 

on  the  other  hand,  usually  employs  As-rich  conditions,  which  should 

result  in  a  large  concentration  of  Ga  vacancies.  For  the  MBE  n-GaAs 

material  used  in  this  study,  the  growth  process  is  only  slightly  As  rich 

[72].  Shown  in  Fig.  4.12  is  the  DLTS  spectrum  of  electron  traps  in  the 

MBE  samples  grown  in  the  first  few  runs  of  the  Riber  system.  The  0.83  eV 

is  not  observable.  Although  there  has  been  a  report  on  electron  traps  in 

MBE  GaAs  [70]  ,  we  believe  that  the  DLTS  signal  in  our  system  has  a  much 

better  signal  to  noise  (S/N)  ratio,  which  allows  better  estimation  of 

2 

energy  levels  as  well  as  concentration.  The  Arrhenius  plot  of  T  t  VS 
1000/T  for  each  trap  is  shown  in  Fig.  4.13.  The  capture  cross  section 
for  each  trap  obtained  from  emission  rate  measurement  is  listed  in  Table 
4.3.  Although  the  electron  trap  of  0.83  eV  is  not  observable  in  LPE  and 
MBE  n-GaAs  material,  we  can  produce  this  tra  . evel  by  changing  the 
stoichiometric  ratio  near  the  surface  layer  as  a  result  of  annealing. 
Shown  in  Fig.  4.14  is  the  electron  trap  spectrum  in  MBE  n-GaAs  sample 
encapsulated  with  SiC^  and  annealed  at  800°C  for  1/2  hr.  As  expected, 
the  trap  MS (0.83  eV)  becomes  observable  due  to  the  annealing  effect  on 
the  encapsulated  sample.  Shown  in  Fig.  4.15  is  the  spectrum  for  LPE 
n-GaAs  sanqple  encapsulated  with  Si^N^  and  annealed  at  800°C  for  1/2  h. 

It  is  surprising  that  a  small  amount  of  the  0.83  eV  trap  level  is  de¬ 
tected.  A  detailed  study  of  annealing  effect  on  MBE  and  LPE  n-GaAs 
material  is  being  undertaken  in  this  laboratory  and  will  not  be  addressed 
here. 

We  believe  that  the  above  experimental  results  support  the 
hypothesis  of  electron  trap  V7  being  related  to  the  Ga  vacancy.  To 


Arrhenius  plots  of  I  t  vs  1000/T  for  electron  traps  in  MBE  n-GaAs  under  various 
experimental  conditions. 
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Table  4.3.  Thermal  Activation  Energies  and  Capture  Cross  Sections 
of  Electron  Traps  in  MBE  n-GaAs. 
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spectrum  of  electron  traps  in  SiO„  encapsulated  MBE  n-GaAs  sample  annealed 
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Fig.  4.15  DLTS  spectrum  of  electron  traps  in  Si  N  encapsulated 


further  test  this  hypothesis  we  have  carried  out  experiments  utilizing 
ion  implantation  as  a  research  tool. 

4.1.4  Study  of  electron  trap  V7  in  ion  implanted  VPE  n-GaAs. 

As  mentioned  in  Section  2.4,  ion  implantation  has  been  a 
powerful  tool  for  semiconductor  research  because  of  its  doping  uniform¬ 
ity,  controllability,  and  reproducibility. 

Shown  in  Fig.  4.16  is  the  spectrum  of  electron  traps  in  VPE 
GaAs  implanted  with  B+  and  annealed  with  Si^N^  encapsulation  at  800°C 

for  1/2  hr.  The  acceleration  voltage  is  250  keV  and  the  dose  is 
14  "2 

10  cm  .  The  projected  range  R  and  its  standard  deviation  AR  can  be 

p  P 

obtained  from  the  tables  in  Appendix  D.  The  peak  concentration  is 

18  -3 

estimated  to  be  2.5x10  cm  .  It  has  been  reported  that  the  carrier 
compensation  due  to  the  implanted  B  atoms  is  not  significant  for  such 
high  dose  implants  and  high  temperature  anneal  [73] .  We  do  not  observe 
significant  conpensation  effect  from  the  I-V  characteristics  shown  in 
Fig.  4.17(a),  or  in  C-V  characteristics  of  the  Schottky  barrier.  There¬ 
fore,  the  net  donor  concentration  is  assumed  to  be  approximately  the 

15  -3 

original  doping  1.2JtlO  cm  .  The  decrease  of  the  trap  concentration 

compared  with  Fig.  4.9  may  be  due  to  Boron  incorporation  into  Ga 

vacancy  sites.  In  order  to  be  sure  that  the  implanted  ions  fit  into  the 

Ga  substitutional  sites,  we  carried  out  Be  implants  into  n-GaAs  [74] 

(Fig.  4.18).  The  implantation  energy  is  250  keV  and  the  dose  is 
14  -2  + 

10  cm  .  The  detailed  planar  p  n  GaAs  diode  fabrication  procedure  has 


been  described  elsewhere  [74] .  Although  the  observed  junction  region 
is  about  1.2  um  [44,  75]  away  from  the  surface,  the  net  donor  concentra¬ 
tion  and  the  trap  concentration  previous  to  the  implant  can  still  be 


Temperature 

spectrum  of  electron  traps  in  B  implanted  \ 


Pig.  4.17  I-V  characteristics  of  Schottky  barriers 

(a)  B  implanted  VPE  n-GaAs; 

(b)  Ga  implanted  VPE  n-GaAs. 
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Fig.  4.18  DLTS  spectrum  of  electron  traps  in  Be  implanted  VPE 


assumed  to  be  about  the  same  as  that  in  the  surface  layer.  Again  a 

significant  decrease  of  the  V7  trap  concentration  is  observed  compared 

with  Fig.  4.9.  As  for  the  trap  V6  (0.58eV),  it  has  been  reported  that 

this  trap  is  related  to  the  radiation  damage  caused  by  Be  implant  into 

VPE  n-GaAs  [76] .  We  have  strong  evidence  from  Fig.  4.9  and  Fig.  4.10 

that  this  trap  does  exist  in  the  unimplanted  material.  Therefore,  the 

defect  responsible  for  V6  can  result  from  the  annealing  of  unimplanted 

material  as  well  as  from  radiation  damage.  We  have  also  performed  Ga 

13  -2 

implants  into  n-GaAs  with  250  keV  implantation  energy  and  2x10  cm 
dose.  As  shown  in  Fig.  4.19,  we  observed  a  very  small  negative  peak 
around  the  temperature  for  the  V7  trap  at  the  present  rate  window 
e  *  47.2  sec  For  a  Schottky  barrier,  only  the  majority  carrier 
(electron)  traps  can  be  observed  if  the  Schottky  barrier  is  good.  Shown 
in  Fig.  4.17(b)  is  the  1-V  characteristics  of  the  Schottky  bacrier  on 
Ga-implanted  n-GaAs.  The  forward  threshold  voltage  is  larger  than  that 
in  Fig.  4.2  and  the  reverse  I-V  characteristic  shows  hysteresis.  This 
larger  threshold  voltage  together  with  the  small  negative  DLTS  peak  may 
be  due  to  the  imperfect  Schottky  barrier  caused  by  the  heavy  ion  implant 
Nevertheless,  the  DLTS  results  on  the  samples  implanted  with  B,  Be  and 
Ga  convince  us  that  the  trap  V7  is  probably  related  to  the  Ga  vacancy. 

As  mentioned  previously,  it  has  been  conjectured  that  level 
V7  is  related  to  oxygen  in  n-GaAs.  There  has  also  been  a  contradictory 
report  that  this  level  is  not  related  to  oxygen  [77] .  This  latter  con¬ 
clusion  is  drawn  from  the  fact  that  the  V7  trap  concentration  by  DLTS 
is  greater  than  the  oxygen  content  estimated  by  SIMS  measurement.  It 
would  be  difficult  to  draw  this  conclusion  if  the  oxygen  content 
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Fig.  4.19  DLTS  spectrum  of  electron  traps  in  Ga  implanted  VPE 
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happened  to  be  larger  than  the  trap  concentration.  To  investigate  this 

controversy,  we  have  done  oxygen  implants  into  GaAs  with  implantation 

14  -2 

energy  of  250  keV  and  dose  of  10  cm  in  order  to  see  whether  the  oxygen 

implant  has  any  effect  on  this  level.  It  has  been  reported  [78]  that 

oxygen  inplants  into  n-GaAs  would  result  in  high  resistivity  material  at 

anneal  temperatures  higher  than  500°C.  There  is  also  a  report  [79]  that 

oxygen  inplants  will  not  necessarily  create  high  resistivity  layers  when 

performed  on  high  quality  epitaxial  GaAs  material.  Indeed,  we  did  not 

observe  significant  semi-insulating  properties  in  the  I-V  characteristics 

shown  in  Fig.  4.20(a)  and  the  C-V  characteristics.  Shown  in  Fig.  4.21  is 

the  DLTS  spectrum  of  the  sample  implanted  with  oxygen  and  annealed  at 

800° C  for  1/2  hr.  Comparing  this  figure  with  Fig.  4.9,  we  observe  that 

the  implanted  oxygen  content  does  not  affect  the  V7  trap  concentration, 

and  we  therefore  conclude  that  oxygen  is  not  involved  in  V7,  the  main 

electron  trap  in  VPE  GaAs.  In  order  to  see  whether  the  incorporation 

of  group  V  elements  affects  this  deep  level,  an  implant  of  N  into  GaAs 

14  -2 

with  250  keV  and  10  cm  was  performed.  The  implanted  sample  was 
encapsulated  with  an  SiO^  layer  and  annealed  at  800°C  for  1/2  hr.  Again, 
we  have  good  Schottky  barriers  with  I-V  characteristics  shown  in  Fig. 
4.20(b)  and  do  not  observe  the  compensation  effect  previously  reported 
[80].  The  DLTS  result  shown  in  Fig.  4.22  indicates  that  the  trap  con¬ 
centration  increases  significantly  compared  with  Fig.  4.9.  This  may  be 
due  to  the  fact  that  nitrogen  fits  into  As  substitutional  sites  as  a 
resonance  state  with  respect  to  the  T  minimum  [81].  Therefore,  the 
effective  concentration  of  the  Ga  vacancy  increases  and  the  V7  trap 
concentration  increases  proportionally. 
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Fig.  4.21  DLTS  spectrum  of  electron  traps  in  0  implanted 
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Fig.  4.22  DLTS  spectrum  of  electron  traps  in  the  N  implanted  VPE  n-GaAs. 


In  the  study  of  the  main  electron  trap  in  VPE  n-GaAs  we  have 


examined  the  annealing  behavior  of  this  deep  level  and  the  effect  of  the 
implanted  species.  In  the  annealing  study,  we  found  that  the  trap  concen¬ 
tration  change  is  consistently  related  to  the  Ga  vacancy  concentration. 

In  the  implantation  study  we  observed  that  the  trap  concentration 
changes  depending  on  whether  the  implanted  species  is  incorporated  into 
Ga  or  As  substitutional  sites.  From  these  DLTS  results,  we  believe  at 
this  time  that  the  main  electron  trap  (V7)  probably  involves  the  Ga 
vacancy. 

4. 2  Studies  of  Electron  Trans  in  .VPE  n~G^./x  Ternary  Alloy 
The  GaAs  P  ternary  alloy  has  been  used  extensively  in  the 

f  X  X 

visible  light  emitting  diode  (LED)  industry.  The  addition  of  nitrogen 

doping  in  this  alloy  yields  a  marked  improvement  as  electroluminescence 

(EL)  efficiency  in  the  indirect  energy  bandgap  region  [82] .  Although  the 

radiative  recombination  mechanism  in  the  GaAs.  P  alloy  or  GaAs.  P  :N 

alloy  system  has  been  extensively  studied  [81-85] ,  there  are  only  limited 

reports  [24-28]  on  the  deep  levels  in  this  material  by  the  deep  level 

transient  spectroscopy  method.  The  nonradiative  recombination  centers 

in  GaAs.  P  alloy  have  been  believed  very  important  to  the  luminescence 
X  —  x  X 

efficiency  [25,26,28,33].  In  this  chapter,  we  will  examine  the  electron 
traps  in  the  GaAs,  P  alloy  system  for  several  compositions. 

X  *X  X 

4.2.1  Material  preparation  and  device  fabrication 

All  the  vapor  phase  epitaxial  (VPE)  GaAs.  P  crystals  used  in 

X  “X  X 

this  work  were  provided  by  the  Monsanto  Company  (St.  Louis)  and  were 

grown  by  an  open  tube  AsH^-pH^  vapor  transport  method  [86].  For  alloy 

+ 

compositions  X  <  0.50  the  layers  were  grown  on  n  (100)  -oriented  GaAs 


substrates,  while  for  X  >  0.50,  n+  (100)  -oriented  Gap  substrates  were 
used.  During  the  first  10  to  50  |j,m  of  growth,  the  ternary  composition 
x  is  graded  at  1%/jsn  from  that  of  the  substrate  to  the  desired  composi¬ 
tion  x  in  order  to  minimize  the  effects  of  lattice  mismatch.  Even  with 
this  composition  grading,  however,  misfit  dislocations  appear  in  the 
epitaxial  layer.  The  indirect-gap  (x  >  0.46)  GaAs1  P  crystals  have 
been  doped  with  nitrogen  during  the  final  10  nm  of  growth  by  adding 

NKLj  to  the  dopant  gas  mixture  [86].  The  N  doping  in  these  crystals  is 

19  -3 

estimated  to  be  ~  10  cm  [86],  The  GaAs,  P  wafers  used  in  this  work 

1-x  x 

were  of  compositions  40+17.,  65+17.,  85+17.  and  1007.  p  on  the  column  V  sub- 

17  -3 

lattice.  The  x  ■  0.40  epilayer  was  doped  with  Te  up  to  ~  10  cm  ,  while 

16  -3 

the  others  were  doped  with  S  up  to  ~  5.5-7.5X10  cm  .  The  "nitrogen 
free"  samples  used  in  this  study  were  prepared  by  etching  away  the  N 
doped  layer  using  a  GaP  etchant  made  by  Transene  Co.  The  etch  rate  is 
~1.5pn/min  at  80°C.  An  etch  period  of~8  min  is  appropriate  for  obtain¬ 
ing  a  virtually  N  free  epilayer,  as  checked  by  photo luminescence  [87], 

All  the  N-doped  samples  were  etched  to  remove  -'■'3-*>  of  the  surface 
layer  to  avoid  composition  variation  [83] . 

Ni  and  Au-Ge  eutectic  alloy  were  evaporated  on  the  backside 

O  O 

of  the  processed  substrate,  with  layer  thickness  of  1000A-2000A,  and 
then  sintered  at  450°C  for  1  minute  to  make  ohmic  contacts.  A1  of 

o  o 

~100QA-2000A  was  evaporated  on  the  epitaxial  GaAs^_xPx  layer  to  form 
Schottky  barriers.  Two  Schottky  barriers  were  cleaved  from  the  same 

sample  and  mounted  on  the  same  To-18  header  used  in  the  DLTS  system 
with  the  two-diode  method  described  in  Section  3.1. 


Electron  traps  in  VPE  n-GaP  have  been  reported  by  several 
authors  [24-28],  and  summarized  by  Lang  [52]  and  by  Mircea  et  al.  [53]. 
Additional  information  is  required,  including  capture  cross  sections, 
along  with  the  emission  rate  to  make  sure  that  deep  levels  observed 
under  different  experimental  conditions  correspond  to  the  same  defect. 

In  this  section,  we  will  report  the  thermal  activation  energies  as  well 
as  the  capture  cross  sections  for  the  traps  observed  under  various 
conditions  in  our  study. 

Shown  in  Fig.  4.23  is  the  DLTS  spectrum  of  electron  traps  in 

N-doped  VPE  n-GaP.  The  thermal  activation  energies  are  as  indicated  in 

Table  4.4.  The  capture  cross  sections  ao  obtained  from  Eq.  (2.13)  are 

* 

also  listed  in  Table  4.4.  Here,  it  is  assumed  that  m  *  0.365  m  and 

e 

mdn  *  1. 20 me  according  to  the  equation  [88]. 

m  “(0.35  +  0.015x)  m  (4. 4) 

and  m^n  *  1.20  mg  for  X  minimum 

The  energy  levels  are  in  good  agreement  with  those  reported  by  Tell 
et  al.  [27],  However,  the  calculated  capture  cross  section  is  much 
larger  than  that  obtained  by  the  capture  rate  measurement  at  the  observed 
temperature  [27].  This  might  be  explained  by  multiphonon  emission  [33] 
(MFE)  according  to  Eq.  (2.12),  in  which  an  activation  energy  for  the 
carrier  capture  is  required.  As  shown  in  Fig.  4.24,  it  is  interesting 
to  notice  that  the  levels  E^  an<*  E^  are  not  detectable  or  resolvable 
in  the  N-free  sample,  and  the  level  E^^  does  not  have  the  same  energy 
as  E^.  We  have  examined  several  samples  and  obtained  the  same  results. 
In  order  to  see  whether  this  change  is  due  to  the  addition  of  nitrogen, 


Temperature  (K) 

Fig.  4.24  DLTS  spectrum  of  electron  traps  in  N-free  GaP. 


106 


18  -3 

N  implants  at  200  keV  with  N  peak  concentrations  of  1x10  cm  and 
18  -3 

5x10  cm  were  performed.  Shown  in  Fig.  4.25  is  the  spectrum  for  a  N 

18  -3 

implanted  sample  with  1x10  cm  peak  concentration  the  energy  level  for 
is  estimated  to  be  0.23  eV,  which  is  between  the  activation  energies 
of  E^  and  If  the  levels  E^,  E^*  and  E^  were  indeed  from  the 

same  origin  but  modified  by  the  existence  of  N  doping,  it  would  be  worth¬ 
while  and  interesting  to  establish  the  relation  between  the  activation 
energy  shift  and  the  N  doping.  The  level  E^  is  estimated  to  be  0.57  eV. 
Levels  similar  to  E^  and  E^  in  the  N-doped  sample  are  not  detected  in 
the  N- implanted  samples.  This  could  be  due  to  the  fact  that  the  N 
concentration  in  the  implanted  sample  is  not  as  high  as  in  the  N-doped 

sample.  Shown  in  Fig.  4.26  is  the  spectrum  for  the  sample  implanted 

18  -3 

with  N  peak  concentration  of  5x10  cm  at  200  keV  and  annealed  at  850° C 
in  Ar  ambient  for  1/2  hr  with  Si02  encapsulation.  The  level  E^  is 
estimated  to  be  0.24  eV.  The  activation  energy  for  E^  is  obtained  as 
0.55  eV.  There  are  two  broad  peaks  which  are  probably  due  to  radiation 
damage.  The  energy  levels  for  these  broad  peaks  have  not  been  deter¬ 
mined.  Similar  broad  peaks  have  been  observed  in  1  MeV  electron  ir¬ 
radiated  Gap  samples  [24],  It  seems  that  there  are  peaks  between  E^ 
and  E^  In  spectrum.  These  peaks  might  be  related  to  the  levels 
E^2  and  E^^  observed  in  the  N-doped  samples.  Further  studies  should  be 

made  in  order  to  verify  that  E  and  E, _  are  related  to  the  existence  of 

12  13 

high  N  doping  instead  of  being  material  dependent. 

4.2.3  Electron  traps  in  VPE  n-GaASq  35 

Shown  in  Fig.  4.27  is  the  DLTS  spectrum  for  N-doped 
GaAs,.  ,.Pn  QC.  The  activation  energies  as  well  as  the  capture  cross 
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Fig.  4.28  DLTS  spectrum  of  electron  traps  in  N-free 


sections  for  the  levels  E._,  Ec/, 

53  54 


and  E^  are  listed  in  Table  4.4.  The 
capture  cross  section  calculation  is  based  on  Eq.  (2.13)  and  Eq.  (4.4). 
Presently  we  cannot  resolve  the  energy  levels  for  E^  and  E,^*  We  can 
estimate  their  activation  energies  by  assuming  the  parameters  in 
Eq.  (2.13)  such  that  [66] 

en  =  k  e‘A®/kBT  (4.5) 

where  k  is  a  constant  for  all  levels  in  this  sample,  i.e.  assuming  the 
same  capture  cross  section  an  .  With  these  assumptions,  the  energy 
levels  for  E^  and  E^^  would  be  0.19  eV  and  0.22  eV,  respectively.  As 
shown  in  Fig.  4.28,  levels  corresponding  to  E  -E^  in  the  N-doped 
samples  are  again  not  observable  in  the  N-free  sample. 

4.2.4  Electron  traps  in  VPE  n-GaAs^  35^0  65 

Shown  in  Fig.  4.29  is  the  DLTS  spectrum  of  electron  traps  in 
N-doped  GaAsQ  ^Pq  Two  overlapping  levels  are  observed,  and  we  can¬ 

not  resolve  their  activation  energies.  However,  as  shown  in  Fig.  4.30, 

the  energy  level  for  EQ1  is  estimated  to  be  0.41  eV  in  the  N-free  sample. 

ol 

Again,  the  spectrum  in  the  N-doped  sample  is  different  from  the  N-free 

sample.  The  evidence  supports  that  E_-  and  EQ1  may  be  from  the  same 

origin.  The  energy  level  of  Eg^  is  very  close  to  that  obtained  for  the 

sulfur-related  deep  level  [28],  which  was  not  detected  in  Gap  but  was 

observable  in  the  alloys  with  composition  0.20  <x<  0.65.  Shown  in 

Fig.  4.31  is  the  spectrum  for  samples  implanted  at  200  keV  with  N  peak 

18  -3 

concentration  of  1x10  cm  and  annealed  at  950°C  for  1/2  hr  wLth  Si02 
encapsulation.  The  levels  EQ1  and  E  0  are  very  similar  to  the  levels 
Ey^  and  E^  in  the  N-doped  samples.  The  level  E^  is  estimated  to  be 
0.26  eV  below  the  conduction  band.  There  is  a  small  peak  of  0.59  eV  in 


Trap 

Activation  Energy 
(eV) 

Capture  Cross- 
Section  (cm2) 

Temperature 

Range  (°K) 

E11 

0.24 

2. 07x10“ 16 

135-160 

E12 

0.27 

3.28xl0‘17 

160-185 

E13 

0.46 

1.49x10" 14 

205-240 

E21 

0.19 

1.95x10" 17 

125-160 

E31 

0.23 

9.35X10"17 

130-160 

E32 

0.57 

-14 

2.52x10 

250-285 

E41 

0.24 

2. 07x10" 16 

135-160 

E42 

0.55 

4.0x10“ 15 

255-295 

E53 

0.25 

3. 24x10" 16 

135-160 

E54 

0.32 

5. 13x10“ 16 

170-200 

E55 

0.52 

3.20x10“ 16 

265-320 

E81 

0.41 

4.38xl0"14 

180-210 

E91 

0.26 

7. 06x10" 16 

135-165 

E93 

0.59 

9.89xl0"14 

245-280 

D1 

0.13 
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the  medium  temperature  range  and  an  even  smaller  peak  at  the  high  tem¬ 
perature  range.  These  two  peaks  are  assumed  to  be  related  to  the  N- 

implant.  Shown  in  Fig.  4.32  is  the  spectrum  for  samples  implanted  at 

18  -3 

200  keV  with  N  peak  concentration  of  5x10  cm  .  We  can  see  that  there 

is  a  very  broad  peak  due  to  radiation-induced  damage.  The  levels  E^ 

and  E^2  are  assumed  to  be  the  same  as  E^  and  E^  in  Fig.  4.31,  because 

they  occur  at  the  same  temperature  for  a  given  rate  window. 

4.2.5  Electron  traps  in  VPE  n-GaAs_  ,  _P_  .. 

- - 0.60  0.40 

Since  GaAs^  gg^O  40  is  a  direct  bandgap  alloy  and  the  epilayer 
is  doped  with  Te  without  the  addition  of  N,  we  expect  a  DLTS  spectrum 
different  from  that  for  the  indirect  bandgap  alloys,  in  terms  of  activ¬ 
ation  energies  and  capture  cross  sections.  As  shown  in  Fig.  4.33  the 
energy  levels  are  estimated  to  be  0.13,  0.18  and  0.29  eV  for  D^-D^.  The 
capture  cross  sections  are  calculated  from  Eq.  (2.13)  by  assuming 


m  “0.089  m  and  m,  *  0.089  m  following  the  equation  [89] 
e  an  e 

m  “(0.068  +  0.052  x)  m 


(4.5) 


and  m  “  m, 
dn 

for  the  T  band  minimum. 

In  this  section  we  have  examined  the  electron  traps  in  VPE 
n-GaAs^_xPx  alloys  with  several  compositions.  It  is  observed  that  the 
DUS  spectra  for  the  unimplanted  GaAs^^P^  alloys  are  similar  to  each 
other.  It  would  be  interesting  to  establish  a  relationship  among  them 
and  a  unified  theoretical  explanation  for  the  variations  in  the 
activation  energies.  There  are  differences  in  DLTS  spectra  among  the 
N-doped,  N-free,  and  N-implanted  samples  for  the  same  alloy  composition. 


The  existence  of  nitrogen  might  play  an  important  role  in  the  formation 


concentration. 
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Fig.  4.33  DLTS  spectrum  of  electron  traps  in  Ga, 


of  deep  levels  in  GaAs1  P  alloys.  Deep  levels  caused  by  the  radiation- 
induced  damage  are  observable  in  the  N-implanted  samples  even  annealed 
at  950°C  for  1/2  hr  with  Si02  encapsulation.  This  study  is  intended 
as  an  initial  survey,  which  is  necessary  for  further  investigations  of 


this  material. 


5.  SUMMARY  AND  CONCLUSION 


In  Chis  work,  detailed  studies  of  the  various  types  of  deep  level 
transient  spectroscopy  systems  have  been  made.  A  special  DLTS  system  has 
been  developed,  employing  a  two-diode  method  to  balance  the  diode  leakage 
currents  at  any  temperature  range.  Using  this  DLTS  system,  we  have 
studied  the  main  electron  trap  in  VPE  n-GaAs.  The  results  suggest  that 
this  deep  level  involves  a  Ga  vacancy.  We  have  also  examined  the  electron 
traps  in  VPE  GaP  and  several  GaAs^_xPx  ternary  alloys.  In  the  following 
sections  the  results  and  conclusions  of  this  work  are  summarized,  and 
suggestions  for  further  investigations  are  noted. 

5.1  DLTS  Studies 

The  typical  DLTS  system  using  a  conmercial  capacitance  meter  (shown 
in  Fig.  3.1)  is  simple  to  operate.  It  is  also  easy  to  calibrate  the 
magnitude  of  the  capacitance  transient.  This  system  has  the  drawback  of 
pulse  distortion  caused  by  the  pulse  transformer  in  addition  to  the  slow 
system  response  time.  The  system  shown  in  Fig.  3.2  eliminates  the  use  of 
the  pulse  transformer.  Pulse  distortion  is  critical  in  the  capture  rate 
measurement  where  a  pulse  with  fast  transition  is  required.  The  slow 
response  time  of  a  commercial  capacitance  meter  together  with  the  system 
overloading  problem  during  the  pulse  period  make  it  necessary  to  gate-off 
the  signal  during  the  initial  period  of  the  capacitance  transient  in 
order  to  avoid  spurious  electrical  effects.  As  a  consequence,  the 
capacitance  transient  after  gating-off  loses  significant  information  if 
the  gate-off  period  is  comparable  to  the  periodic  pulsing  period.  Special 
caution  and  correction  for  data  analysis  are  required  if  a  lock-in 
amplifier  is  used  to  preset  the  rate  window.  In  a  system  employing  a 


dual-channel  boxcar  averager,  it  Is  possible  to  avoid  problems  arising 
from  this  gate-off  requirement  by  setting  the  first  sampling  instant  to 
occur  after  the  gate-off  period. 


In  this  work,  the  gate-off  effect  in  the  lock-in  amplifer  case  has 
been  studied  in  detail.  The  basic  operation  of  a  lock-in  amplifier  for 
the  DLTS  system  has  been  described  in  Appendix  A.  Three  operation  modes 
are  suggested,  dealing  mainly  with  the  lock-in  phase  setting  which  is 
critical  to  the  data  analyses,  including  thermal  emission  time  constant 
and  trap  concentration  calculations.  The  maximum  signal  mode  and  the 
truncated-transient  phase  reference  mode  show  considerable  shifts  in  the 
rate  window  relationship  rmav,/T0  as  a  function  of  lock-in  frequency  as  <" 

9 

shown  in  Fig.  3.14.  The  bias-pulse  phase  reference  mode  is  much  less  •[ 

k 

frequency  sensitive,  and  this  is  clearly  the  best  choice  for  those  who  “ 

► 

prefer  lock-in  operation.  In  a  system  employing  the  dual-channel  boxcar 

averager,  the  measurements  are  typically  made  using  a  relatively  wide  '■ 

,  j 

gate  rather  than  the  generally  assumed  narrow  gate,  because  the  signal-  £ 

to-noise  (S/N)  ratio  is  better  for  larger  gate  widths.  The  thermal 
emission  time  constant  corresponding  to  the  temperature  where  the  DLTS 

* 

peak  occurs  can  be  solved  exactly  from  Eq.  (3.3).  As  a  relatively  good  * 

approximation  to  this  exact  calculation,  the  midpoint  of  the  gate  can  be 
used  to  calculate  the  thermal  emission  time  constants  by  Eq.  (3.4). 

The  improved  DLTS  system  shown  in  Fig.  3.5  was  originated  from  that 
invented  by  Lang  [8],  and  has  a  fast  response  time  of  l-2fisec.  Using 
this  system,  the  gate-off  effect  in  the  lock-in  amplifier  case  is 
negligible  because  the  gate-off  time  is  about  the  same  as  the  pulse 

-j 

duration,  which  is  small  compared  to  the  typical  lock-in  period.  However, 


it  is  difficult  to  use  this  system  to  study  devices  which  have  large 
leakage  currents.  The  leakage  current  is  highly  temperature  sensitive, 
and  is  often  encountered  in  implanted  junctions  if  the  post-implant 
annealing  condition  is  not  adequate.  A  DLTS  system  using  a  two-diode 
method  has  been  developed  for  such  functions.  The  experimental  results 
show  that  the  two-diode  version  is  essential  for  accurate  measurements 
on  diodes  exhibiting  large  leakage  currents. 

5.2  Studies  of  Electron  Traps  in  GaAs 

A  variety  of  electron  traps  in  VPE  n-GaAs  have  been  reported  by 
several  authors  [22,23,52-55].  The  dominant  trap  with  an  activation 
energy  of  0.83eV  has  not  yet  been  identified  conclusively.  This  level 
has  been  conjectured  to  be  related  to  either  a  Ga  vacancy  or  an  oxygen 
impurity.  The  capture  cross  sections  we  obtained  from  emission  rate 
and  capture  rate  measurements  further  support  the  model  proposed  by 
Majerfeld  et  al  [61]  that  the  transition  is  between  the  deep  level  and  the 
L  minima  instead  of  the  T  minimum.  The  trap  concentration  decreases  in 
the  unencapsulated  sample  after  anneal  at  700-750°C  for  \  hr  ,  while  it 
is  unchanged  in  the  unencapsulated  sample  after  anneal  at  750°C  with  As 
over  pressure  for  \  hr.  The  trap  concentration  increases  in  the  Si^N^  or 
SiC^  encapsulated  sample  after  anneal  at  800°C  for  £  hr.  From  isothermal 
and  isochronal  anneals  under  various  ambient  or  encapsulations,  this 
dominant  trap  concentration  seems  to  be  related  to  the  Ga  vacancy. 

Samples  implanted  with  B,  Be,  Ga,  0,  and  N  are  also  investigated. 


From  the  B,  Be,  and  Ga  implants,  we  find  that  the  trap  concentrations  are 
much  lower  in  the  implanted  and  annealed  samples  than  those  in  the  un¬ 
implanted  samples  which  are  annealed  at  the  same  condition.  These  results 


are  apparently  due  to  the  fact  that  the  implanted  species  tend  to  occupy 
Ga  vacancy  sites  substitutional ly.  For  the  0  implant,  we  find  that  the 
trap  concentration  does  not  change  significantly  compared  to  the  unimplant¬ 
ed  sample  through  the  same  annealing  process.  This  indicates  that  this 
dominant  trap  is  not  related  to  0.  We  find  that  the  trap  concentration 
increases  with  N  implantation,  apparently  due  to  the  incorporation  of  N 
into  As  vacancy  sites.  From  this  study,  we  suggest  that  the  dominant 
electron  trap  in  VPE  n-GaAs  is  related  to  the  Ga  vacancy. 

5.3  Studies  of  Electron  Traps  in  GaAs^  Alloys 

Studies  of  electron  traps  in  GaP  have  been  reported  by  several 
authors  [24-28,52,53].  However,  the  discrepancy  of  the  energy  levels 
obtained  under  different  experimental  conditions  is  still  unresolved.  The 
lack  of  experimental  information  in  GaP  as  well  as  the  GaAs1_xpx  ternary 
alloy  motivated  us  to  do  further  studies  by  the  DLTS  method.  For  GaP,  we 
find  that  the  DLTS  spectra  are  different  for  N-doped  end  N-free  samples. 

From  N  implanted  material  we  observe  a  level  with  activation  energy 
between  those  in  the  N-free  and  N-doped  samples.  Further  study  is  necessary 
to  verify  whether  these  levels  are  from  the  same  physico-chemical  origin, 
perhaps  modified  by  the  presence  of  nitrogen.  There  are  two  electron 

traps  observed  in  the  N-doped  samples  but  not  detectable  in  the  N-free 
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sample  or  the  N-implant  sample  with  1 X  10  cm  peak  concentration.  It 

is  possible  that  these  two  traps  are  overwhelmed  by  the  radiation  damage 

in  the  N-implanted  sample.  In  that  case,  these  two  levels  might  be 

related  to  the  N  doping.  We  also  find  deep  levels  related  to  radiation- 

18  .3 

induced  damage  in  the  sample  implanted  at  200  keV  with  5  X  10  cm  N 
peak  concentration,  even  after  950°C  anneal  for  ^  hr.  with  Si02 


encapsulation 


For  the  GaASg  35  alloy,  we  also  find  that  the  DLTS  spectra  for 

the  N-doped  and  N-free  samples  are  different.  For  GaAsrt  ,,P A  ,  c,  a  level 

with  activation  energy  close  to  that  previously  reported  as  a  sulfur- 

related  deep  level  [28]  is  observed  in  the  N-free  sample.  However,  this 

level  is  not  clearly  resolved  in  the  N-doped  GaAs_  ocPft  ,  c,  and  is  not 

detectable  in  the  GaP  and  GaAs^  ^Pq  g^  material  which  is  doped  with  S  as 

a  shallow  donor.  We  also  observe  deep  levels  in  this  alloy  associated 

with  radiation- induced  damage  in  the  sample  implanted  at  200  keV  with  N 

18  -3 

peak  concentration  of  Sx  10  cm  after  high  temperature  anneal.  For  the 
GaASQ  3qPq  40’  *-s  a  direct  bandgap  material,  a  DLTS  spectrum  with 

levels  shallower  than  those  in  Gap  is  obtained.  Although  the  deep  levels 
in  GaAs.  pv  alloys  are  quite  divergent,  there  are  common  features,  such 
as  the  similarity  in  the  DLTS  spectra  for  the  unimplanted  samples  and  the 
difference  between  the  N-doped  and  N-free  samples  with  the  same  alloy 
composition. 

5.4  Suggestions  for  Further  Investigations 

As  pointed  out  in  Section  4.1.3,  the  main  electron  trap  V7  (0.83eV) 
in  VPE  n-GaAs  is  not  detected  in  the  typical  LPE  and  MBE  n-GaAs.  However, 
it  is  observed  in  the  LPE  and  MBE  samples  annealed  at  800°C  for  ^  hr  with 
Si^N^  or  SiC^  encapsulation.  A  detailed  annealing  study  should  be  made  to 
establish  a  unique  explanation  of  this  level  in  epitaxial  n-GaAs.  It  is 
also  interesting  to  see  whether  this  level  will  be  created  by  changing  the 
stoichiometic  ratio  of  Ga  to  As  during  the  MBE  crystal  growth. 

The  role  of  oxygen  in  GaAs  and  GaASj_xPx  is  not  clear  yet.  It  is 
worthwhile  to  do  a  systematic  study  on  oxygen  implanted  samples  with 


available  alloy  compositions  by  the  DLTS  and  photoluminescence  methods  to 
obtain  the  thermal  and  optical  properties  of  oxygen  in  the  alloys.  The 
effect  of  nitrogen  doping  on  the  deep  level  behavior  in  GaAs^  should 
also  be  investigated  further.  It  would  be  a  great  achievement  to 
establish  a  theoretical  explanation  for  the  deep  level  formation  in  this 
compound  alloy,  based  upon  a  combination  of  DLTS  results  with  the  exten¬ 
sive  optical  data  of  Wolford  et  al.  [81,87]. 
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APPENDIX  A 


Lock-in  Amplifier  Operation  and  Numerical  Computation 
Program  for  the  DLTS  System 


The  PAR  HR-8  lock-in  amplifier  has  a  tuned  amplifier  to  pick 
up  the  fundamental  Fourier  component  of  the  input  signal,  a  mixer  which 


generates  a  square  wave  weighting  function  to  multiply  with  the  input  funda¬ 
mental  component,  and  an  output  stage  to  take  the  time  average  of  the  mixer 
output.  There  is  a  phase  control  to  adjust  the  relative  phase  of  the  input 
fundamental  component  to  the  mixer  weighting  function.  As  depicted  in 
Fig.  A. 1,  it  is  easy  to  understand  the  operation  of  a  lock-in  amplifier. 

The  Fourier  components  of  a  periodic  function  f(t)  over  the 

interval  -  4  T  <  t  <  -jj-  T  are 
2  o  —  2  o 


or 


where 


and 


a 

n 


b 

n 


t  f-V2  £(t,COS^  dt’ 
f  f°n 

°  -t  n  0 

o 


c 

n 


.  hza/2 

+  v  > 


n 


=*  tan 


-v!s> 

Vb J 


n 


Figure  A.l  shows  a  square  wave  input,  its  fundamental  component, 
and  the  mixer  outputs  with  various  lock-in  phase  settings.  The  fundamental 
component  is  f.  (t)  *  c.sin  ■  —  sin  .  If  we  multiply  f  (t)  with  a 


Phase  180' 


Phase  270' 


weighting  function  w(t)  in  phase  with  the  input,  the  time  average  of  the 


mixer  output  f^(t)  w(t)  is 

T  /  2 

S  -  “  f  °  f.(t)  w(t)  dt 

°  To  n  1 * * * S 

o 

8 

"  2 
TT 

The  HR-8  lock-in  amplifier  has  a  calibrated  square  wave  to  cali¬ 
brate  the  amplification  gain.  If  the  input  square  wave  has  the  rms  value 
of  its  fundamental  component  equal  to  the  lock-in  sensitivity  setting,  the 
lock-in  output  is  full  scale.  The  peak  to  peak  value  of  a  square  wave  is 
2.22  times  the  rms  value  of  its  fundamental  component.  Therefore,  the  rms 
value  of  f^(t)  is  2/2.22  «  0.901.  The  amplification  gain  G  can  be  calcu¬ 
lated  as 

G  x  /  0.901  -  1  (full  scale) 

TT 

G  =  1.11 


If  the  input  signal  is  a  periodic  exponential  train  as  shown  in 


Fig.  A. 2, 

f(t) 


t  t  +  T  /2 

exp(-  )  exp  ( - ~ — ) 


T 

-f<t< 


T 

o 

2 


The  coefficients  c^  and  0^  are 

-t,  -T 


j  *  T  2  ^  i/2 

exp  (““)  [exp  (-£)  -  Ij/tC^r)  +  ] 

T  T  iT 


tan 


1  T 

-1  o 

2ttt 


If  the  mixer  weighting  function  is  set  as  shown  in  Fig.  A.2,  the  output 

So(td,9,T)  i3 

T  /2 

sn j >  ® *T  )  “  =T-  I*  °  c  sin(|^  +  0  )  w(t)dt 


Before  Gating-Off 


■After  Gating-Off 


2 

n 


exp(-td/T)[l  -exp(-To/T>] 

[  (To/2t>2  +ttV/2 


cos  (0^  +  9). 


(A.  1) 


9 

Where  9  ■  —  x  360 
o 

The  output  signal  S  is  a  function  of  t,,  9,  and  t.  Since  t,  is 

o  a  a 

fixed  and  t  is  a  function  of  trap  energy  level  and  temperature,  only  9  can 
be  adjusted  freely.  There  are  three  convenient  ways  to  adjust  lock-in 
phase  setting — maximum  signal  mode,  truncated  transient  phase  reference 
mode,  and  bias-pulse  phase  reference  mode.  The  first  mode  corresponds  to 
adjusting  the  lock-in  phase  to  give  a  maximum  DLTS  peak  for  each  lock-in 
frequency.  This  mode  can  be  easily  understood  by  setting  9  =  -0^,  i.e.  the 
weighting  function  w(t)  is  in  phase  with  the  fundamental  component.  The 
second  mode  is  obtained  by  setting  w(t)  in  phase  with  the  exponential  train 
after  gate-off,  i.e.  t0  *  0.  In  the  third  mode  one  sets  w(t)  in  phase  with 
the  leading  edge  of  the  bias  pulse  or  the  exponential  train  before  gating 
off,  i.e.  t0  -  -td. 


dS, 


0  for 


max 


The  DLTS  peak  will  occur  at  the  t  value  for  which  — tH 
r  max  dri- 

each  mode.  During  mathematical  manipulation,  one  has  to  notice  that 
cos(0^  +  9)  is  a  function  of  t.  Generally,  it  is  the  numerical  solution  to 


2t 


T 

exp(-To/T){[(Y^  +  2^2t^2  +  Tl2')+  273cos(9^!,1)+IIsin<^04^i^ 


2t  T  T 

-  +  tt2^+  27,]cos(®*P1)-  tts in (9+0^)  -  0 

o 


(A.  2) 


The  appropriate  numerical  computation  program  is  given  in  Section  A. 3. 
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A. 2.  Obtaining  AC  from  Lock-in  Peak  Height 

The  analog  voltage  output  from  the  capacitance  bridge  is 

AC 

r—  exp(-t/r ),  where  S  is  the  conversion  factor  in  pF/Volt.  As  discussed 
SB  B 

in  Section  A. 1,  the  actual  lock-in  peak  output  is  1.11  x  S  (t  ).  AC  can 
*  r  e  o  max 

be  calculated  from 


(k £  N 

1.11  x  S  (t  )  x  (  S,  /  S_  )  -  VT/ 
o  max7  \  B  L J  L 


where  ST  Is  lock-in  sensitivity 

Xi 

V  is  DLTS  peak  height 
L 


then 


SL  *  SB  *  VL 

AC  *  11.1  X  S  (T  ) 
o'  max 


A  simplified  case  is  obtained  for  each  mode  if  t ,/Tq  «  1.  For  the  first 


mode  where  9  *  -0^, 


t  «  0.3485  T 

max  o 

9  =  -0L  «  -24.5° 

S  (t  )  «  0.1738 
o'  max 

and  AC  *  0.518  x  ST  x  S_  x  VT 

L  D  Li 

For  the  second  mode  where  0*0 

T  -  0.424  T 

max  o 

S  (t  )  *  0.1607 
o  max 

and  AC  ■  0.560  x  ST  x  S_  x  VT 

L  D  Lt 

The  third  mode,  where  t_  *  -t.  *  0,  is  the  same  as  the  second  mode. 

vi  a 


A • 3  Numerical  Computation  Program 


00100 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

87388 

m 

01300 

01400 

01500 

01600 

01700 

0l800 

01400 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

O32OO 

O33OO 

03400 

O35OO 

O36OO 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04800 

04900 

05000 


TO  EXCUTE  THIS  PROGRAM  YOU  HAVE  TO  TYPE  DLTSLO.FOR.SYS: IMSL/SEAR 
EXTERNAL  F 

THIS  PROGRAM  IS  TO  CALCULATE  THE  OPTINUM  THERMAL  EMISSION 
TIME  CONSTANT,  PHASE  ANGLE  PHI.  AND  NORMALIZED  LOCK-IN 
OUTPUT  FOR  EACH  OF  THE  THREE  MODES 
MAXIMUM-SIGNAL  MODE  IS  THE  CASE  PHI  =  -THETA, 

TRUNCATED-TRANSIENT  PHASE  REFERENCE  MODE  IS  THE  CASE 
THETAsO . 

glgS  PULSE  PHASE  REFERENCE  MODE  IS  THETA=-DELAY/PERIOD 


$iuiSx  ?  J^t^UM°TIMECCONsff  CAP  :  CAPACITANCE,  PHI 
ANGLE  OF  FIRST  HARMONIC,  THETA  :  PHASE  SETTING 
COMMON  A , BN , XL , XR , EPS , N8lG , XAPP , ITMAX, PI 


:  PHASE 


WRITE(5, 10) 

FORMAT V  DELAY  =  GATE-OFF  PERIOD  =  ' ,$) 

READ(5, 20)  DELAY 
FORMAT (G) 

WRITE(5, 30) 

FORMAT V  LOCK  IN  PERIOD  =  ',$) 

READ (5, 20)  PERIOD 
WRITE (5. 50 ) 

FORMAT ( '  Nsl  FOR  FIRST  MODE,  N=2  FOR  SECOND  MODE 
N=R  FOR  THIRD  MODE,  INPUT  N=  ',$) 

READ (5,20)  BN 
DELAYsDELAY-. 05 
DO  300  1=1,41 

DELAYsDELAY^. 05 
AsDELAY/PERIOD 
XL  =0 . 

XR=3-0 
ITMAX=200 
EPS =.0001 
NSIG=4 

CALL  ZFALSE ( F, EPS , NSIG , XL , XR , XAPP , ITMAX , IER ) 

TYPE  60 , IER 

FORMAT ( f  IER  =  ’,14) 

TYPE  70.  DELAY,  PERIOD,  A 
PRINT  70,  DELAY,  PERI06,  A 

FORMAT ( *  DELAY  =  ',  F7-3/'  PERIOD  =  ',  F8.3/’  A= 
TYPE  80,  XAPP 
PRINT  80,  XAPP 
FORMAT ( '  XAPP  =  '.G) 

TAUMAX  =  PERI0D/(2*XAPP) 

TYPE  90,  TAUMAX 
PRINT  90,  TAUMAX 
FORMAT ('  TAUMAX  =  * .G) 

PHI=ATAN(XAPP/PI)» 180/PI 

CAP  IS  THE  LOCK-IN  OUTPUT  FOR  UNIT  TRANSIENT 
IF (BN -2.0)  150,160,170 
THETAa-1 .*PHI 
GO  TO  200 


=3.0 

«AX=200 


F8.3/’  A=  • ,F6.4) 


THETAsO . 

GO  TO  200 

THETA=-3bO*A 

CAP=2/PI»EXP 


-2  *A*XAPP )• ( EXP (-2  *XAPP)-1 . ) /(XAPP«*2+PI **2 ) 


CAP=2/Pl*EXP (-2 *A*XAPP )•  EXP 
•»0. 5*C0SD(THETA+PHI ) 

TYPE  100,  THETA.  PHI,  CAP 

PRINT  100,  THETA,  PHI,  CAP 

FORMAT (’  THETA  =  ’,G>'  PHI  = 

CONTINUE 

STOP 

END 


' ,G/’  CAP  = 
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B. 1  Dual-channel  Boxcar  Averager  Operation 


In  a  DLTS  system  employing  a  dual  channel  boxcar  averager,  tran¬ 
sient  signals  are  fed  into  the  integrators  with  sampling  instants  at  t^  and 
t2«  For  a  finite  gate  width  w  the  output  for  each  channel  is 

~t,  +  w 


Mo-jr1  . 


S2(t)  = 


=  -  r 

w 

Jt 


t2+"  ,-^dt 


The  output  signal  is 


S (t )  -  S1(t)  -  S2(t) 


-  L(l-e"w/T)(e“tl/T  -e"t2/'r) 


Hence  the  DLTS  peak  will  occur  at  t  where 

max  dT  t 

n 

t  is  the  numerical  solution  to 
max 


(B.l) 


■  0.  Therefore 


e-  (t1-t2)/T  [1+(ti/T  ).e"w/T  (1+w/t  +  t L /T ) ] 

(B.2) 

-[l+(t2/T)-e  W/^T  (1+w/t +t2/T>]  ■  0 

The  numerical  computation  program  is  given  in  Section  B.3. 

B.2  Obtaining  AC  from  the  Boxcar  Peak  Height 

The  PAR  164  integrator  has  a  sensitivity  of  lOOmV  for  a  full 

—  t/T 

scale  output  of  10  volts.  If  the  capacitance  transient  is  AC e  ,  then 
the  DLTS  peak  (volts)  should  be 


x  0.1  -  x  S(T  ) 

10  Sg  max 

where  S_  is  the  conversion  factor  of  the  capacitance  bridge.  AC  can  be 
B 

calculated  from 


For  a  simplified  case  where  w  ■  0, 

Tmax  *  (t2"  V  /  Xn(t2/t1) 

^Tmax^  can  ca^cu^atet^  from  equation  (B.l) 


lumerical  Computation  Program 


i-:* 


00100 

00200 

00300 

00400 

00500 

00600 

00700 

OO8OO 

00900 

01000 

01100 

01200 

01300 

01400 

01500 

01600 

01700 

0l600 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02600 

02900 

03000 

03100 

03200 

O33OO 

03400 

03500 

O36OO 

03700 

03600 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04600 

04900 

05000 

05100 

05200 

O53OO 

05400 

05500 

05600 

05700 

05600 

05900 


C 

c 

c 

c 

c 


30 

40 


50 

C 

C 


60 


70 


80 

90 

100 

150 

200 
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TO  EXECUTE  THIS  PROGRAM  YOU  HAVE  TO  TYPE  EX  DLTSBO.FOR.SYS: 

BOXCAR 

IMSL/SEAR 

THIS  PROGRAM  IS  TO  CALCULATE  OPTIMUM  TIME  CONSTANT,  BOXCAR  DLTS 
PEAK  FOR  FINITE  GATE  WIDTH 
EXTERNAL  F 

COMMON  A . W , T 1 , T2 , XL , XR , EPS , NSIG , XAPP , IIMAX , IER 
WRITE(5,30) 

FORMATS  ENTER  THE  RATIO  A  =  T2/T1  ',$) 

READT5.40)  A 
FORMATiG) 

ITMAX= 100000 
EPS =. 001 
NSIG =4 

DO  200  1=1,10 
P=0.1«I 
TYPE  50.  P 
PRINT  50,  P 

FORMAT ( '  TAU1  NEGLECTS  PULSE  WIDTH.  TAU2  ASSUMES  MIDPOINT, 
TAUMAX  IS  TRUE  VALUE,  GATE  WIDTH  PERCENT  =  ’,G) 

THE  FOLLOWING  CALCULATION  ASSUME  T1  FROM  0.1  TO  50  (ANY  UNIT) 

DO  150  J=1 , 500 

TC  IS  TO  CHECK  T1  AND  T2  ARE  WITHIN  THE  RANGES 
T1=0 . 1*J 
W=P*( A-1 )*T1 
T2=A*T1 

TAU1 =(T2-T1 )/AL0G(T2/T1 ) 

TAU2=(T2-T1 )/ALOG( (T2+W/2. ) /(T1+W/2 . ) ) 

XI  =1  ./TAU1 
X2=1./TAU2 
TYPE  60,  XI,  X2 

FORMAT ( '  XL  SHOULD  BE  LESS  THAN  X2,  XR  GREATER  THAN  XI  ', 

'  XI  =  •  G,'  X2  =  ' ,G/) 

XL=X2/1 .05 

XR=X1*1 .05 

TYPE  70  ,  XL,  XR 

FORMAT ( '  XL  =  ',G,*  XR  =  » ,G) 

TC=(50-W)/A 
IF(TI.GT.TC)  GO  TO  200 

CALL  ZFALSE ( F , EPS , NSIG , XL , XR , XAPP , ITMAX, IER ) 

TYPE  80. XAPP 
PRINT  80,  XAPP 
FORMAT ('  XAPPP  =  • ,G) 

TYPE  90,  IER 
FORMAT (f  IER  =  ' 

TAUMAXsl ./XAPP 
TYPE  100.  T1 ,  TAUMAX.  TAU1 ,TAU2 
PRINT  100,  T1 ,  TAUMAX,  TAU1,  TAU2 

FOMWTj^TI  =  ^  F8.4.3X,'  TAUMAX  =  ',  F7.3,33X,'  TAU1  =',F7.3, 

CONTINUE  "  ’ 

CONTINUE 

END 

FUNCTION  F(X) 

COMMON  A, W,T1 

F=EXP(( A-1 )*T1*X)*( 1+T1*X-EXP(-W*X)*( 1+W*X+T1*X))-(( 1+A*T1* 

X ) -EXP  C-W*X)*( 1 +W*X+A*T 1 *X ) ) 

RETURN 

END 


f, 


,G) 
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00100 

00200 

00300 

00400 

00500 

00600 

00700 

oo6oo 
00900 
01000 
01100 
01200 
01300 
01400 
01500 
01600 
01700 
OlSOO 
01900 
02000 
02100 
02200 
02300 
02400 
02500 
02600 
02700 
02800 
02900 
03000 
03100 
3200 
300 
*00 
>00 
>00 
[00 
300 
03900 
04000 
04100 
04200 
04300 
04400 
04500 
04600 
04700 
04600 
04900 
05000 
05100 
05200 
05300 
05400 
05500 
05600 
05700 
05600 
05900 
06000 
06100 
06200 
oe 
0( 
06500 
06600 
06700 
06&00 
06900 
07000 


1588 


C 

C 

c 

c 


50 

60 

70 

80 

100 

200 

300 

500 


250 

1 

530 

550 

601 

751 

801 

2 

701 

901 

951 


TO  EXECUTE  THIS  PROGRAM - TYPE  'LSFLOC' 

THIS  PROGRAM  USES  THE  ORTHOGONAL  POLYNOMIAL  TO  LEAST  SQUARE 
FIT  THE  DATA  ON  THE  ARRHENIUS  PLOT  AND  CALCULATE  THE  ENERGY 
LEVEL  AS  WELL  AS  THE  CAPTURE  CROSS  SECTION 
COMMON /POL Y /NTERMS ,B(10),C(10),D(10) 

DIMENSION  X(10).F(10).W(10) ,ERROR( 10 ) , PJM1 ( 10 ) , PJ( 10 ) 
DIMENSION  T( 10) ,  FREQ( 10 ) ,TAU( 1 0 ) .ASLOPE ( 10 ) 

REAL  A(80) 

WRITE(5. 50) 

FORMAT ('  ENTER  ONE  LINE  COMMENT  ',$) 

READ (5.60)  A 
FORMAT (80A1) 

PRINT  70 

FORMAT ( *  •»*H*»»m*i»*i*»#»*H*»***t*««»»mm***»') 

PRINT  60. A 
F0RMAT(1X.80A1) 

WRITE  (5,100) 

FORMAT;  '  NPOINTs  NUMBER  OF  THE  DATA  POINTS  s' ,$) 

READ (5, 200)  NPOINT 
FORMAT (G) 

DO  1  1*1, NPOINT 
W(I)*1. 

WRITE (5. 300)  I 

FORMAT; '  T( *  II , * )  =  TEMPERATURE 
READ (5, 200 )T; I) 

X(I)s 1 ,/T(I) 

WRITE (5, 500) I 
FORMAT; *  FRBQC.I1. 

READ (5, 200)  FREQ(I) 

TAU(I)s.424/FREQ(I) 


',$) 


')  =  FREQUENCY  OF  LOCK-IN  ',$) 


TTTAUsT(I)**2*TAU(I) 

XX=1000«X(I) 

PRINT  250 , FREQ(I) ,T(I ) .TTTAU , XX 
FORMAT; '  FREQs  »  ,F6. 1 , f  HZ  TEMPs  ',F8.2,' 
F10.4  •  T**2*SEC  1000/T  s  ' F8.4) 
F(I)4.625E-5»AL0G(t(I)»*2*TAU(1)) 
WRITE(5,530) 

FORMAT; '  THE  EFFECTIVE  MASS  OF  ELE 


KV 


T»«2»TAU 


_  _  ELECTRON  AND  HOLE  IN  SILICON 

C  ARE  1.968,  .5587  5  IN  GAAS  ARE  .068,  .5  '/) 

WRITE(5. 5§0) 

FORMAT; '  EFFECTIVE  MASS  =  *,$) 

RE AD ( 5 , 200 ) EFFMAS 
NTERMS=2 

CALL  ORTPOL (NPOINT.  X,F,W. ERROR, PJM1.PJ) 

TYPE  601,(J,B(J),C(J),D(J  .Jsl, NTERMS) 

FORMAT ( 1  J  s  ',I^('  B(j5=' ,6l6.8,*  C 

C  (J)='  E16.8,'  D(J)=* ,E16.8/) 

TYPE  751 

FORMAT; '  BIISIN(I)  IS  THE  EMISSION  PROPORTIONAL 
C  CONSTANT  WHICH  SHOULD  BE  CONST  ',/) 

PRINT  751 

DO  2  I si, NPOINT 

ASLOB7l)*(pfl)+D(2)«B(1)-D(1))«T(I) 

TYPE  80 1,1, ASLOPE (I ) 

FORMAT (’  ASLOPE II ,’)  =  »,E16.8) 

print  6oi,i,aslop£(i 

CONTINUE 

C0NSTsEXP((D(2)*B( 1 )-D( 1 ))/8. 625E-5) 

TYPE  701 .CONST 
PRINT  701,  CONST 

FORMAT; '  CONSTANTs(3K/M»)*».5*  NC*SIGMAs' ,E16 . 8) 
SIGMAsC0NST»3. 071 966E-22 /EFFMAS 
TYPE  901 .SIGMA 
PRINT  901, SIGMA 

FORMAT; '  CAPTURE  CROSS  SECTION  3  » ,E 16 . 8/) 

TYPE  951, D(2) 

FORMAT ('  THERMAL  EMISSION  ENERGY  3  *,E16.8/) 

PRINT  951 ,D(2) 

END 


SUBROUTINE  ORTPOL ( N PO INT , X , F , W , ERROR , PJM 1 , PJ )U 


COMMON/POLY/NTERMS 
DIMENSION  X(10),W( 
DIMENSION  S(10),T( 
DO  9  J=1,NTERMS 
B( J)=0 . 

D( J )=0 . 

SjjjsO. 


w  ah  a  9  a  «  tr  }R)  uanun  •  t  uri  i  y  t  u  /  u 

RMS, B( 10),C(10),D(10) 
,Wn0),F(l0).ERR0R(10),PJM1(10),PJ(10)U 
,  T ( 1 0 ) , FREQ ( 1 0 ) , TAU (10) 


D&1  ^n&oInt 

ERROR(I)=F(I )-D( 1 ) 

IF(NTERMS  .EQ.  1)  RETURN 
B(  1  )=B( 1 )/S( 1 ) 

DO  12  1=1 , NPOINT 
PJM1 (I)=  1  . 

PJ(I)=X(I)-i}M) 

J=1 
J = J+1 

DO  21  1=1 .NPOINT 
P=PJ(I)»W(I) 

D(J)=D(J)+ERROR(I)*P 

P=P*PJ(I 

Bfj)=B(J)+X(I)«P 
S( J )=S ( J j+P 
D(jj=D(J)/S(J) 

DO  22  1=1, NPOINT 
ERROR(I)=ERROR(I)-D( J)*PJ(I) 

IF(J  .EQ.  NTERMS)  RETURN 
B(J)=B(J)/S(J) 

C(J)aS(J)/S(J-1) 

DO  27  1=1, NPOINT 
P=PJ(I) 

PJ(I)=(X(I)-B(J))*PJ(I)-C(J)»PJM1(I) 
PJM1 (I)=P 
GO  TO  20 


3(J)=BCJ)+X(I) 
3(jj=S(J)+P 
)(jj=D(J)/S(J) 
)0  22  1  =  1  ,NPOI 


-D(J)*PJ(I) 

RETURN 
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TO  EXECUTE  THIS  PROGRAM - TYPE  'LSFBOX* 

THIS  PROGRAM  USES  THE  ORTHOGONAL  POLYNOMIAL  TO  LEAST  SQUARE 
FIT  THE  DATA  ON  THE  ARRHENIUS  PLOT  AND  CALCULATE  THE  ENERGY 
LEVEL  AS  WELL  AS  THE  CAPTURE  CROSS  SECTION 
DIMENSION  X(10),F(10),T(10), TAU (10), ASLOPE (10) 

REAL  A(80) 

WRITE (5, 50) 

FORMAT ('  ENTER  ONE  LINE  COMMENT  ',$) 

READ! 5. 60) A 
FORMAT ( 80A 1 ) 

PRINT  1150 
FORMAT ( ' 

PRINT  1200. A 
FORMAT! IX, 80AI) 

WRITE (5.  IPO) 

FORMAT! '  ENTER  PULSE  WIDTH  IN  usee  =  ',$) 

READ! 5, 1400)  PD 
FORMAT (G) 

PRINT  I5OO. PD 

FORMAT!’  PULSE  WIDTH  =  »,F8.3,»  usee') 

WRITE  (5,100) 

FORMAT!  '  N POINTs  NUMBER  OF  THE  DATA  POINTS  s' ,$) 

READ! 5. 200)  NPOINT 
FORMAT !G) 

DO  1  1=1 .NPOINT 
WRITE!5,300)  I 


1500. PD 

!'  PULSE  WIDTH 


*,F8.3,'  usee') 


FORMAT! 


,11,')  =  TEMPERATURE 


)RMAt('  EN 
SADT5.6OO) 
IITE($.510 


THE  SECOND  SAMPLING  INSTANT  IN  msec 


x  viuiA  axif  /  —  x  x  w  an  x  uii  u  j  y/ 

form1?!o)0)t  11 
X(I)=1./T(I) 

WRITE (5 1 5002 

FORMAT ! '  ENTER  THE  FIRST  SAMPLING  INSTANT  IN  msec  '  ,$) 

READ(5,600)T1 

FORMAT (G) 

WRITE(5, 1000) 

FORMAT! *  ENTER  THE  SECOND  SAMPLING  INSTANT  IN  msec  ',$) 

READi5.600)T2 

WRITE($.510) 

FORMAT! f  ENTER  THE  APERTURRE  DURATION  IN  usee  ',$) 

READ! 5 .520)  AD 
FORMAT! G) 

TllsTI 

T22=T2 

T1sT1-.001*PD+.0005*AD 
T2=T 2- . 00 1 »PD+. 0005*AD 
TAU(I)s.001»(T2-T1)/ALOG(T2/T1) 

TTTAU  =T (I ) **2*TAU ( I ) 

AXs1000*X(I) 

RATIOsEXP . 00 1 «T 1 /TAU ( I ) ) -EXP ! - . 00 1 *T2/TAU (I)) 

PRINT  1 100, Til ,T22,AD,T(I) , TTTAU, AX 

FORMAT!  •  T1=  ',F6.3.'  T2=  ',F6.3,'  MSEC  ','ADs  *  ,F8. 3. 

»  USEC','  TEMPi  *  .P8. 3,  *  K  , '  T”2»TAU=  ',F10.3,?  T**2- 
SEC  ',»  1000/T  s  J ,F8. 4) 

PRINT  1250,  RATIO 

FORMAT!'  AMPLITUDE  RATIO=  EXP(-TI/TAU)-EXP(-T2/TAU)=  ',F7.5) 
FTl)[g8  *  625E-5*AL0G(T(I)**2*TAU(I) ) 

FORMAT? j5THE  EFFECTIVE  MASS  OF  ELECTRON  AND  HOLE  IN  SILICON 
ARE  1.963,  .55587  ;  IN  GAAS  ARE  .068,  .5'/) 

WRITE (5. 556) 

FORMAT!'  EFFECTIVE  MASS  =  ' ,$) 

READ (5. 560)  EFFMAS 
FORMAT (G) 

CALL  LISQUA ( N  PO INT , X , F , SLO  PE , GONST ) 

TYPE  751 

FORMAT!'  EMISIN(I)  IS  THE  EMISSION  PROPORTIONAL 
CONSTANT  WHICH  SHOULD  BE  CONST  ',/) 

PRINT  751 

ASLOPE ( 1 5 If? J I J-GONST ) «X ! I ) 
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TYPE  80 1,1, ASLOPE (I) 

FORMAT ('  ASLOPE? ' ,11 , ’ )  =  ’.E16.8) 

PRINT  801, I, ASLOPE (I) 

CONTINUE  „  , 

CONST =EXP ( -GONST/8 . 625E-5 ) 

TYPE  701, CONST 

FORMAT (•  CONSTANT = ( 3K/M* ) ** . 5*  NC«SIGMAa’ ,E16.8) 


PRINT  701. CONST 
C0N$T«3.C 
01, SIGMA 


SIGMA: 


071986E-22/EFFMAS 


TYPE  99I, SIGMA  „ 

FORMAT (’CAPTURE  CROSS  SECTION  a  ’,E16.8/) 

PRINT  901 .SIGMA 
TYPE  951, SLOPE 

FORMAT ( '  THERMAL  EMISSION  ENERGY  ?  ’,E16.8/) 

PRINT  951, SLOPE 
END 

SUBROUTINE  LISQUA (NPOINT, X,F , SLOPE , GONST ) 

DIMENSION  X(10 ) ,F( 10  ) 

XBARsO . 

YBARaO. 

PRODUCaO. 

SQUAREaO . 

DO  10  lal, NPOINT 
XBARsXBAR+X ( I) 

YBARaYBAR+F(I) 

PRODUC  aPRODUC+X ( I ) • F( I ) 

SQUAREaSQUARE+X ( I ) **2 
CONTINUE 

XBARsXBAR/NPOINT 

YBAR-YBAR/NPOINT 

SLOPEa( PRODUC -NP0INT*XBAR*YBAR)/(SQUARE-NP0INT*XBAR**2 ) 

GONST aYBAR-SLOPE*XBAR 

END 
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APPENDIX  D 


Projected  Range  Statistics  of  Implanted 
B,  Be,  Ga,  0,  and  N  in  GaAs,  Obtained  from  LSS  Theory 
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